Technische Universitat Miinchen 

Physik Department 
Institut fur Theoretische Physik T30d 



The Cosmological Constant and Discrete 

Space-Times 



Dipl.-Phys. Univ. Florian Bauer 

Vollstandiger Abdruck der von der Fakultat fiir Physik der Technischen Universitat Miinchen 
zur Erlangung des akademischen Grades eines 

Doktors der Naturwissenschaften (Dr. rer. nat.) 
genehmigten Dissertation. 

Vorsitzender: Univ.-Prof. Dr. Lothar Oberauer 

Priifer der Dissertation: 1. Prof. Dr. Manfred Lindner, 

Ruprecht-Karls-Universitat Heidelberg 

2. Hon. -Prof. Dr. Wolfgang Hillebrandt 



Die Dissertation wurde am 31.07.2006 bei der Technischen Universitat Miinchen einge- 
reicht und durch die Fakultat fiir Physik am 23.08.2006 angenommen. 



2 



Abstract 

In this thesis the cosmo logical constant is investigated from two points of view. First, we 
study the influence of a time- dependent cosmological constant on the late-time expansion 
of the universe. Thereby, we consider several combinations of scaling laws motivated by 
renormalisation group running and different choices for the interpretation of the renormal- 
isation scale. Apart from well known solutions like de Sitter final states we also observe 
the appearance of future singularities. As the second topic we explore vacuum energy in 
the context of discrete extra dimensions, and we calculate the Casimir energy density as 
a contribution to the cosmological constant. The results are applied in a deconstruction 
scenario, where we propose a method to determine the zero-point energy of quantum fields 
in four dimensions. In a related way we find a lower bound on the size of a discrete grav- 
itational extra dimension, and finally we discuss the graviton and fermion mass spectra 
in a scenario, where the extra dimensions form a discrete curved disk. 



Zusammenfassung 

In dieser Arbeit wird die kosmologische Konstante aus zwei verschiedenen Blickwinkeln 
untersucht. Als erstes behandeln wir den EinfluB einer zeitabhangigen kosmologischen 
Konstante auf die Entwicklung des Universums zu spaten Zeiten. Dabei betrachten wir 
mehrere Skalengesetze, die vom Renormierungsgruppenlaufen herruhren, und auBerdem 
verschiedene Moglichkeiten, die Renormierungsskala festzulegen. Neben bekannten Lo- 
sungen wie dem de Sitter Kosmos beobachten wir auch das Auftreten von Singularitaten 
in endlicher Zukunft. Der zweite Schwerpunkt dieser Arbeit stellt Vakuumenergie in 
diskreten extra Dimensionen dar, wo die Casimirenergiedichte als Beitrag zur kosmolo- 
gischen Konstante berechnet wird. Im Rahmen von Deconstruction verwenden wir die 
Ergebnisse, um eine Moglichkeit zu finden, die Nullpunktsenergie von Quantenfeldern in 
vier Dimensionen zu bestimmen. Ebenso leiten wir eine untere Schranke fur die Grofie 
einer diskreten gravitativen extra Dimension her. Und schlieBlich diskutieren wir die 
Massenspektren von Gravitonen und Fermionen in einem Modell, bei dem die extra Di- 
mensionen eine diskrete gekrummte Scheibe bilden. 
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1 Introduction 



Several years after the discovery of the accelerated cosmological expansion, the question 
what drives this behaviour is still open. Since general relativity with well known matter 
sources like dust and radiation always exhibits a decelerating universe, we have to expect 
something new that explains the observed behaviour, which has gained great support by 
recent observations 121 El EJ E] • The name "dark energy" has become common to 
describe all energy forms that are able to yield the acceleration, although other origins 
are often put into the same category The amount of possible frameworks and models that 
have been proposed to describe DE has grown quite huge as can be seen in recent reviews 
about this topic [Tj El El [TU1 [TT] . Interestingly, the first dark energy candidate was already 
introduced by Albert Einstein almost a century ago. He considered a positive cosmologi- 
cal constant (CC) in order to realise a static cosmos, where the CC compensates the dust 
matter. However, it turned out that this space-time was unstable and additionally not 
consistent with the afterwards observed cosmological expansion. These days, Einstein's 
CC has come back as the simplest candidate for dark energy since it is just a constant 
in the action for classical general relativity. Furthermore, in the context of quantum field 
theory the vacuum or zero-point energy of quantum fields contributes also to the CC. 
But this quantum origin also involves the so-called cosmological constant problem 112; ■ 
which is generally considered to be one of the biggest mysteries in physics. The core 
of this problem is the fact that in quantum field theory the absolute value of the CC 
has not been determined yet, because the corresponding calculation leads to an infinite 
value. Furthermore, naive estimations using energy cutoffs are many orders of magnitude 
above its measured value. Due to this extreme discrepancy and the fact that dark en- 
ergy has become dominant very lately in the cosmological evolution, other explanations 
have been proposed. Apart from new energy forms like, e.g., scalar field condensates 
(quintessence) |T3|, that contribute to the energy content of the universe, modifications of 
the theory of gravity or even of the space-time structure have become subject of intensive 
investigation. Instead of discussing all these possibilities, we will in this thesis concen- 
trate on two subjects in the context of vacuum energy as a major source of dark energy. 
First, we will study a CC that becomes time-dependent due to quantum effects. In this 
framework we discuss the late time cosmological evolution and the occurrence of future 
singularities. The second topic of investigation are space-times with discrete compact 
extra dimensions, where the appearance of zero-point energy and the properties of fields 
will be discussed. 

To introduce the first subject let us start on the classical level, where the CC is just a 
constant term in the action for general relativity, implying that it remains constant in Ein- 
stein's equations, too. Assume for the moment that the current cosmological acceleration 
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is just due to a positive CC in a universe with only cold dark matter (ACDM cosmology), 
then one finds that our space-time is approaching an empty and static de Sitter universe. 
However, other dark energy sources might change the cosmological fate significantly in 
comparison to the ACDM case. For instance, they could cause a final domination of mat- 
ter or even future singularities [H], where space-time collapses to a big crunch or gets 
torn apart in a big rip event ^Hj. Especially the big rip singularity would require some 
very unusual energy forms like scalar fields with negative kinetic terms and other "exotic" 
things. But as we will show in this work, such extreme final states might also appear 
with the CC once it is allowed to be time-dependent or more generally scale-dependent. 
Constants that become scale-dependent occur quite generally in quantum theories like 
quantum electrodynamics, where the fine-structure constant becomes energy- dependent 
by renormalisation group effects. The corresponding renormalisation group equations 
(RGE) describe the "running" of the constant as a function of the renormalisation scale \x 
and masses of quantum fields that are involved in the theory. Similarly, one obtains the 
RGEs for the CC and also for Newton's constant from the effective action of quantum 
fields on curved space-times jTH H3 [TBI QSj. For free fields the RGEs can be calculated 
exactly thereby leading to a scale dependence of the CC as we will see later. Also certain 
theories of quantum gravity [20J |2H 122] can be a source for a scale-dependent CC or New- 
ton's constant. Unlike running coupling constants emerging from interacting quantum 
fields, where the renormalisation scale can be easily interpreted as external momentum or 
temperature, the identification of the scale for the CC is not always given by the theory 
since there is no external momentum. To study the RG running in a cosmological context 
we will therefore explore as candidates for the renormalisation scale several scales that are 
characteristic for the cosmological evolution. This will be the Hubble scale and the sizes of 
the particle and event horizons. In addition, the consequences of several RGEs following 
from different frameworks and motivations will be discussed as well. For each combination 
of RGEs and renormalisation scale identifications we will finally solve Einstein's equations 
and determine the possible final states of the universe. Apart from the pure curiosity of 
the researcher, the results of this analysis might be valuable also for deciding which of the 
above combinations are reasonable after one has accepted or rejected the idea of future 
singularities. 

The second major part of this thesis deals with fact that vacuum energy is also sensi- 
tive to external conditions, which might come from a non-trivial space-time structure. 
The Casimir effect [2B] is a famous example, where the quantum field zero-point energy 
depends on the boundary conditions that are imposed onto the quantum fields. In the 
original setup the boundary conditions follow from two parallel conducting plates in four 
dimensions, where the Casimir effect implies a force on the plates. While in this case 
the corresponding vacuum energy usually cannot be identified directly with contributions 
to the CC, a space-time with compact extra dimensions (ED) might lead to an effective 
four dimensional (4D) CC that depends on the properties of the EDs [24J. Generally, 
one can say that the resulting Casimir energy of massless fields scales inversely with the 
size of the EDs. Here, one encounters a severe problem because small EDs produce large 
contributions to the effective 4D CC in contrast to its observed tiny value. Too large 
EDs, on the other hand, could be easily observed by fifth force experiments and other 
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methods [23 EfE EZ| • One has therefore to cope with stringent bounds on the size of the 
EDs. Instead of working with large EDs one can obtain small CC contributions also by 
considering massive quantum fields, which yield a suppression of the Casimir energy. This 
property allows to keep both the ED and the Casimir energy small. In this thesis we will 
investigate the Casimir effect and the behaviour of quantum fields more deeply in the 
context of discretised EDs. In contrast to continuous EDs the discrete structure implies 
some interesting features like a finite number of Kaluza-Klein modes or typical lattice 
theory attributes. More motivation comes from the fact that a discretised space-time 
structure exhibits a minimal physical length scale that could in principle regularise space- 
time singularities in general relativity and respectively ultra-violet (UV) divergences in 
quantum theories. It therefore represents a useful concept for quantum gravity |2%] . By 
the way, as EDs have not been observed yet, the property of being discretised represents 
a reasonable possibility. One should just imagine some kind of higher dimensional solid 
state physics. Furthermore, it was recently found that discretised EDs can be described 
within a fully 4D context called "deconstruction" jSnHSDj- This model-building approach 
avoids problems that emerge when extra- dimensional field momenta reach the fundamen- 
tal Planck scale, which can be considerably below the 4D Planck scale. Coming back to 
the CC problem, the vacuum energy of 4D quantum fields in a deconstruction setup is, as 
expected, divergent, which just leads to the renormalisation group running effects men- 
tioned above. Thus the overall CC scale is undetermined and has to fixed by observations. 
However, in this work we propose a well defined prescription to assign a finite value to 
the CC by employing the correspondence between deconstruction and discretised higher 
dimensions. We will demonstrate this idea within a specific deconstruction model, where 
also the suppression of the Casimir energy by massive fields will be applied. 

Going one step further by discretising gravity in the EDs one approaches new effects 
in the form of strongly interacting massive 4D gravitons [HU IH2j- This strong coupling 
behaviour leads to an upper limit for masses and energy scales appearing in the theory. 
Interestingly, the limit might depend also on the size of the ED. As a consequence, the 
ability to suppress the Casimir energies by large field masses is reduced since the masses 
have to lie below the strong coupling scale. In this context we will derive some limits 
on the ED size by requiring that the resulting CC does not exceed its observed value. 
Finally, we consider on a more formal level a six dimensional (6D) model, which involves 
a discretised curved disk. In this scenario we discuss the implementation of fermions and 
gravitons, and furthermore investigate the effect of the curvature and its influence on the 
mass spectra of the resulting 4D fields. 

The structure of this thesis is given as follows: in Chap. |2]we introduce a number of scal- 
ing laws for the CC together with some renormalisation scale identifications. Afterwards 
we explain how to solve Einstein's equations with a time-dependent CC and Newton's 
constant and subsequently discuss cosmological late-time solutions and the corresponding 
fates of the universe. Chap. El is devoted the Casimir effect in continuous higher dimen- 
sions, where the vacuum energy for a five dimensional (5D) setup is derived. In Chap. 0] 
we introduce the discretisation of the fifth dimension and explicitly show how to calculate 
the Casimir energy in the corresponding scenario. Several effects like bulk masses and 
lattice artefacts are discussed, too. As an application we apply in Chap. El our results to 
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a deconstruction model and propose a way to determine the zero-point energy of the 4D 
fields. Discretised gravitational EDs are the main subject of Chap. where we derive 
some bounds on the size of the ED. In Chap. [7|we investigate a 6D scenario with two 
curved and discretised EDs. We calculate the corresponding mass spectra of gravitons 
and fermions and close the chapter with a little application. Finally, we summarise the 
results of this work and present our conclusions. 

Conventions 

Throughout this work we use the Einstein sum convention, and we set the speed of light c, 
Planck's constant H and Boltzmann's constant fcs to unity. Other conventions are given 
in the text, where necessary. 



2 Time-dependent Cosmological 
Constant 



2.1 Quantum Effects and the Cosmological Constant 



In this chapter we will investigate a positive CC as the most prominent candidate for 
dark energy. As a component in Einstein's equations of classical general relativity it can 
be treated as a perfect fluid with a constant energy density A > and an equation of 
state uj = p/A = — 1 that corresponds to a negative pressure p = —A. On the quantum 
level the CC emerges as vacuum energy of quantum fields with the same equation of state 
as the classical CC. It is therefore an unavoidable constituent of the matter content of the 
universe. Unfortunately, one does not know how to calculate its value in a unique way, 
because it can be written in the form of a quartically divergent momentum (p) integral 
like J d 3 p ■ p. Respectively, for compact dimensions one obtains the infinite sum of zero- 
point energies \p. In Minkowski space one usually eliminates the infinite vacuum energy 
by normal ordering in QFT since it has no influence on flat space-time physics. Gravity, 
on the other hand, is sensitive to all forms of energy and matter, and we thus have to deal 
with vacuum energy in cosmology. The naive assumption of an UV cutoff to regularise 
the infinite integral at some known energy scale M x leads to an unobserved high value of 
the CC, which illustrates the old CC problem ^2]. For example, let M x ~ 10 3 GeV be 
the energy scale, where supersymmetry is assumed to be broken. This leads via 



to a mismatch of about 60 orders of magnitude. Without extreme fine-tuning in the theory 
the naive cutoff method for determining the CC obviously does not work and should be 
rejected. Fortunately, the procedure of renormalisation in QFT can handle infinities, 
thereby leading to a dependence of the renormalised constants on some energy scale /i. In 
many cases, this renormalisation scale can be identified with an external momentum, or at 
least with some characteristic scale (e.g., the temperature) of the environment. Studying 
QFT on curved space-time fOJ EI! leads to infinities in the effective action or in the 
vacuum expectation values (VEV) of the energy-momentum tensors of the fields. This 
can be treated by renormalisation to yield a scale-dependent or running CC and a running 
Newton constant. However, the absolute values are still not calculable, but the change 
with respect to the renormalisation scale can be calculated via RGEs. In the following 
we will investigate the influence of RGEs originating in QFT jTTJ |33l |33j and quantum 
gravity [20]. Unlike the running coupling constants in the standard model of particles, 
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here, the physical meaning of the scale is not given by the theory since there is no external 
momentum. In the cosmological context that we consider, it is reasonable to identify the 
renormalisation scale with some characteristic scales in cosmology, which will be done in 
Sec. 12.31 For phenomenological reasons, we should request that the scale \x does not change 
too much over cosmological time scales. Once we have fixed the combination of RGE and 
scale identification we are able to discuss in Sec. 12.51 the late-time behaviour and the fate 
of the universe with a running CC. In addition to a running CC we will also consider a 
running of Newton's constant for one case in Sec. 12.61 Please note, that in contrast to 
dark energy scenarios with a time- dependent equation of state as in Refs. here, the 
equation of state of the CC is still exactly uj = — 1. Nevertheless, it is possible to obtain 
an effective time-dependent equation of state [HE] due to a non-standard scaling of the 
matter energy density with the cosmological time. This non-standard scaling is explained 
in Sec. I2.4[ where we discuss how to solve Einstein's equations on a Robertson- Walker 
background when the CC and Newton's constant depend on the cosmological time. 



2.2 Scaling Laws 



According to QFT on curved space-times the CC and Newton's constant G are subject to 
renormalisation group running like the running of the fine-structure constant in quantum 
electrodynamics. The corresponding scaling laws of these constants depend crucially on 
the considered quantum fields and their masses. In the following we consider two different 
RGEs emerging in this framework and in addition one scaling law that emerges in the 
context of "quantum Einstein gravity" .2Q.. Throughout the text we will identify the CC 
with the corresponding vacuum energy density A. 

Let us first start with non-interacting quantum fields on a curved space-time, namely 
a Friedmann- Robert son- Walker universe with a positive CC. For one fermionic and one 
bosonic degree of freedom with masses nip and ms, respectively, the 1-loop effective action 
can be written in the form UJ 



S, 



cff 




4 4 
m F — /Tig 

32tt 2 



Ric 

167TG 

Ric 



A + [ D + In 



16vr 2 



- i 



m F / B \ 
H ) 



C 



(2.2) 



where D = ^Euier + lim n ^ 4 (n — 4) _1 is a divergent term, which does not depend on 
the renormalisation scale \x. Furthermore, £ is a coupling constant 1 , and the variable C 
represents all further terms in the effective action, that are neither proportional to the 
Ricci scalar Ric nor to the vacuum energy density 



A :-- 



A 



8ttG 



1 In the action S — J d 4 x^^g [ 



Ric-2A 



167rG -r 2 ( t ) \a^ a — 2 i m + fR* ]^ J of a scalar field <j) on a curved space- 
time, the constant £ occurs in the coupling term £ • Ric • <fi 2 between the scalar field and the Ricci 
scalar Ric. 



2.2. Scaling Laws 
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The relevant /3-functions in the MS-scheme for the vacuum energy density A and Newton's 
constant G are obtained by the requirement that the effective action S e g must not depend 
on the renormalisation scale /z, 

dS^fj 
fi—— = 0. 
d/x 

Because of this condition, A and G have to be treated as /i-dependent functions in 
Eq. ()2.2)1 . which consequently have to obey the RGEs given by 

dA mt — ran d / 1 \ 1 r / >. in9 1 9i 

Note, that the divergent term D has dropped out, leaving over just the masses ttif/b 
and £. Assuming constant masses, the RGEs can be integrated, hence, the equation for 
the vacuum energy density reads 

A(/i) = A (l-giln^), A := A(// ), (2.3) 

where Ao denotes the vacuum energy density today, when the renormalisation scale fi has 
the value /io- Moreover, the sign of the parameter 

depends on whether bosons or fermions dominate. In this context, a real scalar field 
counts as one bosonic degree of freedom, and a Dirac field as four fermionic ones. The 
generalisation to more than one quantum field in the RGE, can be achieved by summing 
over the fourth powers of their masses. For Newton's constant G we obtain the RGE in 
the integrated form 

G(fi) = °° , G :=<?(/«,). (2-5) 
l-<7 2 ln£ 

Again, we omit the generalisation to more fields, that follows from summing over the 
squared masses of the fields. For one bosonic and one fermionic degree of freedom the 
mass parameter g 2 is given by 

q 2 := ^ [(£ - \)ml - >|] . (2.6) 

Finally, we remark, that Eq. ()2.3j) for the running vacuum energy density A(/i) was derived 
in a renormalisation scheme, which is usually associated with the high energy regime. It 
is therefore not known to what extend it can be applied at late times in cosmology. In 
addition, to avoid conflicts with observations, the field content has to be fine-tuned to 
obtain \qi\ < 0(1), which we assume for the rest of this work. Unfortunately, the corre- 
sponding covariantly derived equations for the low energy sector are not known yet (36j . 
Therefore, we prefer to work with the above RGEs, which were derived in a covariant way, 
and study the consequences and the constraints on the mass parameters q\ and q 2 - 

The second scaling law follows from a RGE that shows a decoupling behaviour. Consid- 
ering only the most dominant terms at low energy, the corresponding /5-function for A is 
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given by /i^p = Aifi 2 , where A\ ~ ±M 2 is set by the masses M and the spins of the fields. 
RGEs of this kind have been studied extensively in the literature, see Refs. [Tl\. Assuming 
constant masses and [Aq to be of the order of today's Hubble scale H one obtains 



A(/i) /i 2 M 

— — = L + -Li — , -Li ~ ± — 
Ao A*o m p 



L ^±T72i (2-7) 



where L := 1 — L\ and Mp denotes the Planck mass today. Here, the running of A is 
suppressed since \L%\ <C 1 for sub-Planckian masses M. 

The last scaling laws come from the RGEs in quantum Einstein gravity [201 EI]- in 
this framework the effective gravitational action becomes dependent on a renormalisation 
scale /i, which leads to RGEs for A and G. An interesting feature is the occurrence of an 
UV fixed-point [37J in the renormalisation group flow of the dimensionless quantities 2 A/x 
and Gfi 2 at very early times in cosmology corresponding to \i — > oo. Motivated by strong 
infrared (IR) effects in quantum gravity one has proposed that there might also exist an 
IR fixed-point in quantum Einstein gravity [38] leading to significant changes in cosmology 
at late times 3 , where /J — ► 0. If this were true, one would obtain the RGEs 

A _ /i 4 G _ $ 

A ~~ /, 4 ' C ~ ;i2' ^ > 

A /i Ct /i 

where fiQ corresponds to A = A and G = Go. In the epoch between the UV and IR 
fixed-points, A and G vary very slowly with /x, and we will treat them as constants in this 
region. Therefore, we assume that the scaling laws ()2.8|) are valid from today until the 
end of the universe. 



2.3 Renormalisation Scales 

In order to study the effects on the cosmological expansion due to the scaling laws of 
Sec. 12.21 we have to define the physical meaning of the renormalisation scale /i. Un- 
fortunately, the theories underlying the scaling laws often do not determine the scale 
explicitly 4 , apart from the usual interpretations as an (IR) cutoff or a scale characterising 
the physical environment (temperature, external momenta). In our cosmological setting 
we will investigate three different choices for /j, given by the Hubble scale H, the inverse 
radius R -1 of the cosmological event horizon, and finally the inverse radius T _1 of the 
particle horizon. 

Let us consider our universe on sufficiently large scales, where it is well described by the 
Robertson- Walker metric given by 

d^d^-a^^ + rW), (2.9) 

2 In this work A denotes the vacuum energy density corresponding to a CC, whereas in many articles 

about quantum gravity A is the CC A itself. 
3 However, it was argued in Ref. that the region of strong IR effects can never be reached. 
4 In the framework of Refs. Newton's constant was found to depend explicitly on the Hubble scale. 



2.3. Renormalisation Scales 
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where k fixes the spatial curvature. In accordance with recent observation we will ignore 
the spatial curvature and consider in the following the spatially flat (k = 0) metric 

ds 2 = dt 2 - a 2 (t)dx 2 , (2.10) 

where the scale factor a(t) depends only on the cosmo logical time t and not on the spatial 
coordinates x. 

Our first candidate for the renormalisation scale // is the Hubble scale 

*( t ):=|>, (2,1) 

which describes the actual expansion rate of the universe. On the other hand, the horizon 
scales R and T describe the cosmological evolution of the future and the past, respectively. 

In universes with late-time acceleration like the ACDM model there usually exists a cos- 
mological event horizon. Its proper radius R corresponds to the proper distance that a 
(light) signal can travel when it is emitted by a comoving observer at the time t: 

R ^--^fwy (2 - 12) 

In the case that the universe comes to an end within finite time the upper limit of the 
integral has to be replaced by this time. Similar to the event horizon of a black hole, 
the cosmological event horizon exhibits thermodynamical properties like the emission of 
radiation with the Gibbons-Hawking temperature [H]. 



The counterpart of R is the particle horizon radius T, which is given by the proper distance 
that a signal has travelled since the beginning of the world (t = 0): 

T(t) -- =a{t) lwy (213) 

In a simple cosmological model, where the universe begins at t — and then evolves 
according to the decelerating scale factor a(t) oc t n with < n < 1, the particle horizon 
radius reads T(t) = t/{l — n). Now we assume that at some time t a de Sitter phase sets 
in, which corresponds to a scale factor a(t) oc exp(H (t — to)) for t > to. This leads to a 
radius function T that grows exponentially with t at late times: 

T{t > to) = \ T(t ) + ±-J exp(H (t - t )) - (2.14) 

Since the asymptotic behaviour of the scales H, R and T plays a major role in this work, 
it is plotted in Fig. 12.11 for a universe with dust-like matter and A and G being positive 
and constant. Note that if one of the horizon radii diverges, one says that the horizon 
does not exist. The geometric meaning of the scales given here has been discussed, e.g., 
in Ref. jl2j, and some arguments for the event horizon in the context of the CC can be 
found in Refs. 
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Figure 2.1: In a spatiafly flat ACDM universe with constant A and G t'he radii of the event 
horizon R and the Hubble horizon H~ l approach the same constant value for t — > 
oo. At early times the particle horizon radius T grows linearly with time and 
exponentially with time for t — ► oo. At the time to = 0, 99 Hj7 (today) the relative 
vacuum energy is given by Qao = 0, 73. 



2.4 General Relativity with Time-dependent Constants 

In this section, we derive the evolution equation for the cosmic scale factor a(t) in the 
framework of the spatially isotropic and homogeneous Friedmann- Robert son- Walker uni- 
verse with a time-dependent CC and Newton's constant. On this background, radiation 
and pressureless matter (dust) can both be described by a perfect fluid with energy den- 
sity p and pressure p = u>p, where the constant u characterises the equation of state. For 
instance, dust-like matter like CDM has the equation of state uj = and incoherent radi- 
ation lu — g, respectively. The corresponding energy- momentum tensor for these energy 
forms reads 

T af3 = (p + p)u a vP - pg af3 \ 

with u a being the four-velocity vector field of the fluid. With our choices of the renormal- 
isation scale, G and A depend only on the cosmic time t. From Einstein's equations 

G a(3 = 8nG(Ag a P + T aP ) 

and from the contracted Bianchi identities G a/3 .p = for the Einstein tensor G a/3 , we 
obtain the generalised conservation equations 

[GAg a/3 + GT a % = 0, 

whose a = component reads 



G(A + p) + G(A + p + 3|p(l + oj)) = 0. 



2.4. General Relativity with Time-dependent Constants 
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Here the dot denotes the derivative with respect to the cosmological time t, and we do 
not assume T.* = 0. For constant A and G the last equation can be integrated to yield 
the usual scaling law for the matter energy density p oc a~ 2<y ® +l \ where we have for our 
convenience introduced the equation of state parameter 5 



Q:=i(l + 3u;). 



(2.15) 



Note that for non-constant A and G this simple scaling rule for the matter content is 
not valid anymore, because it is now possible to transfer energy between the matter 
and the vacuum, in addition to G ^ 0. At this stage we have to admit that this energy 
transfer implies an effective interaction between the gravitational sector (A, G) and matter, 
which is not part of the original Lagrangian. In this sense it should be compared with 
gravitational particle production ^3] resulting also from the interplay of gravity with 
quantum physics. 

Since we have at this point no information about p as function of time or the scale 
factor, we have to combine the Friedmann equations for the Hubble scale H := - and the 
acceleration -, 



f) +4 = ^G(t)(A(t)+ P (t)), 



a 
a 



57T 



G(t)(A(t) - Qp(t)), 



(2.16) 
(2.17) 



in order to eliminate the matter energy density p. The left-hand side of the resulting is 
abbreviated by F(t): 



F{t) :=- + Q 
a 

By introducing the constant 

K :- 



+ 



k 



STY 



GA • (1 + Q). 



(2.18) 



Ho 2 



87rG Ao(Q + 1) n A0 (Q + l) 



(2.19) 



with the relative vacuum energy density ^ao = 8tiGoA /(3Hq) and the Hubble scale H 
at the time t = to, we obtain the main equation in the compact form 



K F(t) 



AG 
A G*o 



(2.20) 



Here, we can now insert the RGEs for A and G from Sec. 12.21 together with a choice of the 
renormalisation scale /i from Sec. 12.31 For the scale identifications p = Rr 1 and p = T _1 
we have to solve integro-differential equations, whereas p = H just leads to an ordinary 
differential equation. 



5 For a dominant matter energy density /)> A and flat spatial curvature, the acceleration quantity q 
is given by the negative value of the Q. 
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2.5 Cosmological Evolution at Late Times 



In this section we study the late-time evolution of the universe with variable A and G, 
thereby assuming from today on the validity of the scaling laws ()2.3J) - (|2.8J1 and the correct 
identification of the renormalisation scale /i with the scales (I2.11j) - (j2.13j) . Our aim is to 
determine in all nine cases the possible final states of the universe. This depends, of 
course, on the choice of parameters, but we restrict ourselves to parameter values which 
comply vaguely to current observations. Today, at the cosmological time to = 0, 99//^ 1 
we fix the initial values H and Q^o = 0, 73 by observations. Furthermore, the initial 
value T of the particle horizon radius would be fixed if the past cosmological evolution 
was known from the Big Bang on. In contrast to this, the value Rq of today's event 
horizon radius depends on the future cosmological evolution and is treated here as a free 
parameter. In the following we derive some properties of the solutions analytically, in 
particular, we study the stability of (asymptotic) de Sitter solutions and the occurrence 
of future singularities, where the scale factor or one of its derivatives diverge within finite 
time US]- For simplicity we denote a big rip or a big crunch by the lowest order 
divergent derivative that is positive or negative, respectively. Since some combinations of 
scaling laws and renormalisation scales lead to complicated equations we derive in some 
cases only approximate or numerical statements. 

In order to solve Eq. (|2.20j) numerically, we have to remove the integrals in the definitions 
(|2.12j) . ()2.13|) of R and T by a differentiation with respect to t. Using the relations 



R = RH-1, T = TH + 1 



an ordinary differential equation for the scale factor can be obtained. Let us, for instance, 
consider fi = R~ l in combination with the RGEs ()2.3)1 and (|2.5J) . First, we solve the main 
equation (l2~20l) for R, 



n{t) 

which has the time derivative 



R(t) 

R 



exp 



K F(t) 



qi - q 2 K F(t) 



(2.21) 



a (gg - qi)K F 
a ( Ql - q 2 K Fy 





I K F-1 1 


■ exp 


_q 1 - q 2 K F_ 



1 

Rq 



(2.22) 



where we have substituted R by Eq. (|2.21J) . This equation can be integrated numerically. 
Afterwards we have to check whether the functions R and T, calculated from the numerical 
solution of a(t), agree with R and T that follow directly from Eq. (|2.20p . If they do not 
match, the numerical solution has to be discarded. Some of these cases are illustrated and 
discussed in Sec. 12.61 Furthermore, solutions involving a negative matter energy density 
are questionable on physical grounds. This happens when the vacuum energy density A 
becomes greater than the critical energy density p c = 3H 2 / (8irG), which follows from the 
first Friedmann equation (|2.16|) . 
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Figure 2.2: Lambert's LF-function is the inverse function of z = xe x . For —1/e < z < it has 
two real-valued brunches an d W-i(z), and for z > only Wo(z) is real. 



2.5.1 A = A (l - qi In £), p = H 

For the given scaling law and scale choice Eq. ()2.20j) reads 

K (H +(Q + l)H 2 ) = 1 - q x In J^. (2.23) 

-no 

We first look for asymptotic de Sitter solutions by applying H = and H — > i/ e . Thus 
the final Hubble scale i? e is given by 

^ = q -V A0 W u ( -^-e'/A. (2.24) 



H Q 2 \<?i^ao / 

Here, W u (z) with it = 0, — 1 denotes one of the two real-valued branches of Lambert's 
W-function, which is the solution of z = xe x , see Fig. 12.21 For qi > there is always one 
solution for H e , and for qi < two solutions exist if the argument of W u in Eq. (|2.24J) is 
greater than — e -1 . This means either q\ < n\ or q\ > n 2 with 

2 2 

These solution are stable if 



d(K H) 



dH 



= -^L-2H c K (Q + l)<0, 



where we used Eq. (|2.23|) . For positive q% this condition is always fulfilled, whereas 
for negative q% it means W u ( -^-^ e 2 ' qi ) < —1, which follows from Eq. ()2.24j) . Again, 
the argument of W u is constrained yielding q% > n 2 - Therefore only q\ > n 2 , which 
includes q% > 0, leads to stable de Sitter solutions. Using the phase space relation (|2.23J) . 
which is plotted in Fig. 12. 3[ we conclude that for other values of q\ the cosmological 
evolution will always end within finite time in a Big Crunch singularity, where H — > 
and H — ► — oo. 
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Figure 2.3: The phase space relation (j'2.23j) is shown for q\ = — o; —3, where the universe ends in 

a big crunch singularity, and respectively for q\ = — T which leads to a de Sitter 

final state {H = 0). In all cases the evolution begins at H = Hq, where Q\q = 0, 73 

and Hq < 0. 



2.5.2 A = A (L + fi = H 



M6 



Using this choice for A and \i the main equation (|2.20J) becomes 



- + H 2 (Q - Li(Q + 1)Oy ) = H 2 L n A0 (Q + 1), 



which has the exact solution 
a(t) 



a 



with the parameter 



sinh([Q + l]{tt A0 L (l - L&^Hot) 
n:= [(g + l)(l-L 1 fi A0 )]- 1 . 



At early times, H t <C 1, it describes a power-law expansion a(t) oc t n , whereas at late 
times it approaches the de Sitter expansion law a(t) oc exp(H c t) with the final Hubble 
scale H e given by 



H e = Hq 



1 — Lif^o 

Further aspects of this case have been studied, e.g., in Refs. [T7j . 



2.5.3 AG = A G 4, /' = H 

Following Sec. l2.2l we assume that this scaling law is valid from today on. Here, Eq. ()2.20j) 
is given by - + BH 2 = with B := Q — (Q + l)fiAo- It has an exact power-law solution 

^ = [(1 + B)(t - h = const., 

a 



that exhibits a constant acceleration da/ a 2 = —B > as well as Q A = const. Note that 
a similar solution was found in Refs. [3%] . 
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§ 0,6 




Figure 2.4: The lower bound i?omin of the initial event horizon radius Rq from Eq. (|2.26j) for 

negative q\. For Rq < i?omin a final de Sitter state does not exist. 



2.5.4 A = A (l - qi In fi = R 



-l 



At first we look for de Sitter solutions with a(t) oc exp(H c t), where the inverse Hubble 
scale and the event horizon radius approach R — > R e = H^ 1 in addition to H = 0. 
Plugging this asymptotic form for a(t) into Eq. ([2.20)1 . one arrives at 

— = q 3 x =l + qi\nx, 

where the variables x := R e / Ro and q 3 := K (l + Q)/Rq > have been introduced. The 
solutions for x are given by 



x 



Re 
Rq 



2q 3 



(2.25) 



involving Lambert's IT-function W u (z) with u = 0,-1 (Fig. 12.2)1 . Since the IT-function 
is real only for arguments z > — e _1 and 2 < for W_i(z) we obtain for negative q\ the 
constraint q 3 < — yexp(— ^- — 1), which implies an lower bound for the initial value Rq 
of the horizon radius as shown in Fig. 12.41 



-Ro > -Romin : = 



— n — 772 ex P 
5liiA0-"n 



2 



(2.26) 



If i?o is smaller than this minimal value, then Eq. (|2.25j) has no positive solutions and 
a final de Sitter state does not exist. For R = -Romin there is exactly one solution x = 
exp(— — i), for higher values Rq there are two solutions. In the case of a positive value 



'/1 
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of the parameter qi, the initial value Rq must be smaller than 1/ a/ H^Qaq. Otherwise the 
final horizon radius R e is smaller than the initial one, x < 1. Both cases are plotted in 

Fig, ma 

Since we have found several de Sitter solutions, we now have to study the stability of 
these final states. Therefore, we write K F as a function of K F = K (H + (Q + 1)H 2 ), 



K F = gi 



1 f l-K F 
H - — exp 



where we used R = RH — 1. In the final de Sitter state we have K F = and R = R c 
1/H e = const. Near this point we can neglect H in the function F and replace H by 



KnF 



MQ + iy 



For a stable solution it is required that 



d(K F) 



d(K F) 



as 



Z Kq 



exp 



1 - K F 

Qi 



< 0, 



(2.27) 



in the final point, where K F = K (Q + 1)H% = q%x 2 . With q 3 and x from above, this 
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yields the stability condition 



2 ^3 2/ 9l 



W H ( — e 



-i 

< 



implying that there are no stable de Sitter solutions for positive values of q%, because 
the ^-function is positive. Thus a big crunch singularity (F, H — > — oo, R — > 0) might 
happen for certain values of initial conditions and parameters, or there exists no solution. 
A big rip singularity (F, H — > oo, R — > 0) does not exist because it is not compatible 
with i2 — > oo. 

For negative gi we get the condition 

W u [^e 2/qi ] > -1, 



9i 

which means that only the solution with u = is stable. This renders the final event 
horizon radius R e = RqX unique. 

Finally, we take a closer look at the ratio R e /Ro as a function of the mass parameter q\. 
For initial values Rq < l/(-£^oV^X\o), which means q% > 1, there is a certain range of 
values of q\ where no solutions for R e exist. This range is again given by the requirement 
that the argument of the Pi^-function must be greater than or equal to — e^ 1 , leading to 
the conditions 

* £ w^xy 01 * a wZpry (2 ' 28) 

In Fig. 12.61 the exclusion range for q\ is obvious for g 3 > 1. In the case that q\ lies above 
this range, the unstable solution for R e is reached first during the future cosmic evolution. 
For qi below this range, the stable solution is nearer to the initial value Rq than the 
unstable one, however, both solutions for R e lie below i? - Initial values R > 1/ (i/ V^Xvo) 
(i.e. 53 < 1) lead to stable final states with R e > Rq for all negative values of q\. Moreover, 
all stable de Sitter solutions are realised except the last one corresponding to the parameter 
region 

2 

i^Ao<l, ?i> Wo[ _ 1/{e i q3) y 

Numerical solutions show that in this last case the universe will end in a big rip singularity, 
where F, H — > 00 and R — ► 0. A big crunch (F — > —00) is not possible for negative q\ 
because it contradicts R — > 0. This can also be seen from the time-derivative of the 
function K F given by 

K F = q 1 !l = q 1 (H-^, q x < 0. 

This becomes positive for H < (big crunch), thus preventing a further decrease of K F 
and the final collapse of the scale factor. 

More on the numerical solutions of this case can be found in Sec. 12.61 and Ref. [33], where 
additionally the running of Newton's constant is taken into account. 
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Figure 2.6: This figure shows the ratio R c /Rq from Eq. (|2.25|) between the final and the ini- 
tial event horizon radius as a function of q\. We plotted four cases for different 
values of q% = 1/(HqRqQaq), where the choice = 0,8; 0,9 corresponds to an 
initial radius Rq > l/(ffoV^Ao)> and q% = 1,1; 1,2 to Ro < 1/(#oV^ao)- I n the 
latter case there are no solutions for a certain range of values of q\ as described by 
Eq. (|2,28|) , The thick lines show the {u = 0)-branch of R e /Ro, and the thin lines 
the (u = — l)-branch, respectively. 



2.5.5 A = A (L + L^), ^ = R' 1 

In the de Sitter limit, H = and R —>■ R c = H~ l , we find here the solution 

R l = (H 2 n A Q - Lii?o)/L , 
but it is unstable for all possible values of L\, which follows from 

= R c HILqVLxo > 0. 

dS 

For L\ < our numerical calculations did not yield any solution that is compatible with 
Eqs. (I2.12j) and ()2.20|) . This is also true in a certain parameter range for L\ > 0, where 
otherwise a Big Rip singularity occurs, where F, H —>■ oo and R — > 0, see Fig. 12.71 A big 
crunch is not possible since K F is bounded from below. 



d(K F) 
d(K F) 
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(a) (b) 

J? 2 

Figure 2.7: Corresponding to Sec. l2.5.51 we plot (a) the function KqF = Lq+L\^ for L\ = ±0, 3 
and (b) the Hubble scale H(t) for L\ = +0, 1; Rq = 1, 1; Oao = 0, 73 and Q = \ 
(dust), where a big rip singularity occurs. 



2.5.6 AG = A Go4, A* = R' 1 

Here, one cannot find a prediction for R e in the de Sitter limit (H = 0, R — > i? c 
i?" 1 ) since Eq. ()2.20|) only leads to the constraint -Rq = 1/(HqQ\ ). To find solutions 
corresponding to other values of Rq we impose the ansatz a(t) = ao(t — £i) n with ao, n., £i = 
const. The Hubble scale is given by H = n^—ti)^ 1 which implies due to H > that t± <t 
for n > and £i > t for n < 0, respectively. For < n < 1 the event horizon does not 
exist, in the other cases R can be calculated exactly via Eq. (|2.12j) : 

R(t) - a(t) r dt ' -(t tr fe-^) (1 " W) -(^^) (1 " n) 

Note that for negative n there is a big rip singularity at t — > ti, where a(t) diverges and the 
universe ends: t^, — t±. For n > 1 the scale factor a(t) describes power-law acceleration 
and there is no future singularity which means £g — > oo. As a result, both cases yield 



n — 1 

Thus Eq. f!2.20|) can be solved exactly, which determines the constant n, see Fig. |2I 



(z - |) ± J(x - i) 2 - x(x - (Q + 1)) 
= V(Q + 1) < 2 ' 29 > 

withx := = R 2 Hfa A0 (Q + l). (2.30) 
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Figure 2.8: The exponents ni2(x) from Eq. (|z.29|) . which occur in the power-law expansion 
scale factor a(t) = ao(t — t\) n , are shown. In Sec. 12.5.51 we find n = n^Ry/ Kq), 
whereas n = ^(Tq/Kq) in Sec. 12.5.51 

For Q > we find (Q + < n-i < 1 for all positive values of x, therefore 712 can be 
dropped as the event horizon does not exist. In the case R$ < 1/(HqQ^ ) the constant n = 
n\ is negative leading to a big rip singularity at the time t — t\, and > 1/(HqQ^ ) 
implies a positive n corresponding to power-law acceleration. 



2.5.7 A = A (l - qi lnf o ), fi = T _1 

Solving Eq. ()2.20|) for this choice of A and \i can be done approximately. First we propose 
the ansatz 

a(t) = a exp[u 2 (t - t^ 2 ] (2.31) 
for the scale factor with the constants a , u and t\. Therefore 

H = 2u 2 {t-t x ), - = Au 4 (t-t 1 ) 2 + 2u 2 , F = 2u 2 + Au\Q + l)(t-t 1 ) 2 , 

(X 

and the particle horizon radius at late times reads 

T = exp(V (i - tx) 2 ) / exp(-u 2 (t' - ti) 2 )dt', t > t x = const. 



exp(w 2 (t — 1\ 



2u 



erf(w(t-ti)) 



J i. 



For t — > oo the series expansion of the square brackets reads C + 0(t ) with C > and 

Zo 

T 



thus T w C exp(w 2 (t— ti) 2 ). In this limit the given ansatz for a{t) solves KqF = 1 — qi In 2a 



exactly and determines the constant u 2 = |gii?o^Ao- Therefore we have found that 
Eq. (|2.31|) represents an approximate late-time solution in the case q\ > 0, which we call 
super-exponential accelerated expansion. 
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Figure 2.9: Here, the numerical and the approximate solutions for the scale factor a(t) and the 
Hubble scale H{t) from Sec. I2.5.7l are shown for negative q\ = —0,2. Furthermore, 
we have used S^ao = 0, 73; TqHq = 3 and Q = ^ (dust). In this case a big crunch 
singularity will occur in the future. Note that the analytical approximation from 
Eq. (|2.32() works very well except near the future singularity. 



The numerical solutions in the case qi < exhibits a future big crunch singularity. How- 
ever, at times sufficiently before this event the scale factor can be well approximated by 
the ansatz 



which implies F 



a(t) = a exp[— u 2 (t — t\ 
-2u 2 + 4u 4 (Q + l)(t-t 1 ) 2 and 



(2.32) 



T = exp(— u 2 {t — t 1 



IT 



—erf(iu(t - ti)) 
Liu 



const. 



We expand the square bracket term around t — > t±, where the scale factor is maximal, 
and find T « Cexp(-u 2 (t - h) 2 ) with C > 0. Solving K F = 1 - qtln f leads to u 2 = 
— jqiHQfl\Q in accordance with qi < 0. The numerical and approximate solutions are 
shown in Fig. | 



2.5.8 A = A (L + M = T' 1 

In the late-time de Sitter limit, a(t) oc exp(H c t), the particle horizon radius has the 
asymptotic form 

T oc exp(H e t) — > oo for t — > oo, 

which follows from Eq. (j2.14|) . This means that K F —>■ L and thus determines H 2 = 
HqQ^oLq. To verify the stability of this solution we calculate 

T T 2 T 2 / 1 

K oF = _2L x i • ± = -2L X ^- (h+± 
i T 2 T 2 \ 7^ 
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At the initial time (today) KqF is negative for L\ > and KqF keeps on approaching the 
de Sitter limit unless a sign change in KqF would occur. But this requires a negative H, 
which is only possible if K F = H + (Q + 1)H 2 < at some time. Since K F > L > 
we conclude that the de Sitter limit will be reached always. For L\ < the initial values 
of KqF and KqF are positive and KqF approaches the de Sitter limit. Again, a sign 
change in KqF is only possible for H < 0, which requires KqF < at some time. These 
conditions cannot be realised because KqF > as long as KqF > 0. In summary all 
reasonable values of L\ lead to a stable de Sitter final state. 



2.5.9 AG = A G 



Ml 
Mo 



2 » 



By using at late times the power-law ansatz a(t) = a Q (t — ti) n with a,Q,t,n = const., we 
find 

F =(n 2 (Q + l)-n)(t-t l )- 2 , 



and the particle horizon radius is given by 



T = a {t-t 1 ) n 



t0 dt' 

L./o ^0 



* dt' 
to W)\ 



t - u 



n 



(t 



1 — n 



(2.33) 



where To = T(tn) > depends on the past evolution of the scale factor. In the case n < 
a positive Hubble scale H = n(t — ti) -1 requires t, to < ti, which implies T < and thus 
the non-existence of the particle horizon. For < n < 1 the second term in Eq. ()2.33|) 
can be neglected at late times, and Eq. ()2.20|) can be solved exactly, leading to 



with x 



(x - \) ± ^(x-i) 2 -x(x-(Q + l)) 



(Q + l) 



T 2 



-^ = T 2 H 2 n A Q(Q + l). 

An 



This is the same equation as in Sec. 12.5.61 with Rq replaced by T , see also Fig. 
Here, solution n 2 is the physical one because of (Q + 1) _1 < n 2 < 1, that corresponds 
to a decelerating universe for t — > oo. In comparison with Sec. 12.5.61 the identification of 
the particle horizon radius as renormalisation scale yields a complementary cosmological 
behaviour. 



2.6 Time-dependent Cosmological and Newton's 
Constant 

In the previous sections we have discussed only the running of the CC except for the 
scaling laws (|2.8|) . where the change of Newton's constant is comparable with that of the 
CC. For the scaling law (|2.3J) we have required that the mass parameter qi is small in order 
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to obtain a viable phenomenology. Since the largest known field masses are of the order 
of 10 2 GeV we had to assume fine-tuning or some unknown suppression mechanism to 
achieve the smallness of q\. What happens to Newton's constant when it is controlled by 
the RGE (|2.5jl ? As we can see from Eq. (|2.6j) . the corresponding mass parameter q 2 is of 
the order m 2 Go = m 2 /Mpj and therefore strongly suppressed by the Planck mass Mpi even 
for larger masses m. Obviously, the running of G is suppressed from the beginning, which 
agrees with the strong bounds on the time- variation of G, see, e.g., Ref. |T7]. Additionally, 
this property has the advantage, that today we are far away from the Landau pole of G(n), 
where the function F diverges in Eq. ()2.2U|) . Since the RGE for G follows directly from 
the effective action ()2.2)1 . we are interested in its influence on the cosmological evolution. 

As an example we will discuss the combination of the RGEs (j2.3|) for the CC and (|2.5|) 
for G that were both derived from the effective action ()2.2|) . As renormalisation scale \x 



we use, according to Eq. (j2.12j) . the inverse of the cosmological event horizon, /i = R 
With these preliminaries Eq. ()2.20|) reads 



l + giln| 
1 + <?2 In Ro 



KoF = — -4-, (2.34) 



leading to the ordinary differential equation ()2.22|) for the scale factor a(t). Unfortu- 
nately, finding explicit solutions of this equation seems to be rather difficult because of 
the strongly non-linear form of the equation. Therefore, we solve it numerically to show 
the characteristic future cosmic evolution corresponding to four different cases, which re- 
sult from the parameter choices q\ = ±2 and q 2 = ±0, 1. These values for q\ mean that the 
relevant mass scale m should be near A^ A ~ 10~ 3 eV. Actually, the only known particles 
with such a low mass are neutrinos. This indicates that the influence of higher mass fields 
is suppressed, or these fields have decoupled, respectively. Note that due to the suppres- 
sion by the Planck scale, the realistic value of q 2 should be much lower than ±0, 1. Here, 
we used a large value for q 2 just to explore the differences due to the sign of q 2 . Concerning 
the differential equation, we fix the initial conditions by using observational results, i.e., 
the relative vacuum energy density is given by fl^ — 0, 73, no spatial curvature = 0) 
and only dust (with an equation of state parameter Q = 0,5) and the CC as relevant 
energy forms in the present-day universe. The acceleration parameter go — ^l(^o) is de- 
termined by Eq. (j2.18|) : q = Q^(l + Q) +Q(Qk — 1)- Today's value of the horizon radius Ro 
is unknown, so we have to estimate it. Since it should be the largest physical length scale 
and the universe seems to be almost de Sitter-like, we assume the horizon radius to be a 
bit larger than the inverse Hubble scale, Rq > Hq . Finally, all dimensionful quantities 
are expressed in terms of today's Hubble scale H (Hubble units). Figures 12 . 1 0H2. 13l show 
the numerical results for different values of the initial radius Rq of the event horizon. The 
graphs in each of the four figures illustrate the scale factor a(t), the Hubble scale H(t), 
the acceleration q(t), the event horizon radius R(t), and F(t) as functions of the cosmic 
time t, respectively. The last graph displays Eq. ()2.34|) . i.e., K F as a function of the 
radius R/Rq. 

The first observation from the numerical solutions is, that for a positive value of qi the 
cosmic age decreases with respect to the age t of the standard ACDM universe, whereas 
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for a negative qi the age increases. In the latter case, we observe only big rip solutions 
and de Sitter final states. However, for positive q\ a big crunch may also occur, but no 
stable de Sitter states exist. For positive values of qi and g 2 (see Fig. 12.13]) . we have not 
observed any big rip solutions. Hence, the final state may be either a big crunch or a forever 
expanding universe, where the Hubble scale approaches a finite positive value, but the 
event horizon radius R goes to infinity. This is a contradiction, because an asymptotically 
constant Hubble scale H > implies a finite event horizon radius R ~ H -1 in the far 
future, which is not the case here. Obviously, this numerical solution is not a solution of 
the original equation ()2.20|) . For qi > and q 2 < (see Fig. l2.12|) the big rip events in the 
numerical solutions occur at a finite and large value of the horizon radius R. Again, this 
behaviour is not compatible with the vanishing of the horizon radius at such an event. 
Therefore, we can reject these numerical solutions, too. For positive qi we thus observe 
that if a solution exists it has to be a big crunch. 

In summary, we conclude that the cosmological final states occurring in the numerical 
solutions in this section do not differ from the case with constant G, discussed in Sec. 12.5.41 
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Figure 2.10: The cosmological evolution for the parameter choice q± = —2 and q2 = — 0, 1 and 
different values of the initial event horizon radius Rq. The fate of this universe is 
either a stable de Sitter state when choosing Rq = 1,20; 1,30, or a big rip (BR) 
in the case Rq = 1, 10; 1, 15. KqF is bounded from above. Nomenclature: Scale 
factor a, Hubble scale H = -, acceleration q = event horizon radius R and its 
: i : ' : ; i value Rq. For the function KqF see Eqs. (|2.H4|) . and (|2,19|h for the mass 
parameters qi,q2 see Eqs. (|2.4|) and ()2.6I) . 
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Figure 2.11: For the parameter choice q\ = —2 and qi = +0, 1 the cosmic evolution is not very 

different from the case q2 = — 0, 1 (Fig. I2,ll7|) , In the future, there is either a stable 

de Sitter state for Rq = 1, 20; 1, 30, or a big rip (BR) when Rq = 1, 10; 1, 15. KqF is 

bounded from below. Nomenclature: Scale factor a, Hubble scale H = -, accelera- 

tion q = |§, event horizon radius R and its initial value Rq. For the function KqF 

see Eqs. I|2.H4|I and (|2.19|) . for the mass parameters qi, qi see Eqs. (|2.4|) and (|2.Hjl . 
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Figure 2.12: The cosmological evolution for different values of today's horizon radius Rq for the 
case q± = +2 and q2 = —0,1. The solutions for Rq = 1,09; 1,10 exhibit a big 
crunch (BC), whereas the initial conditions Rq = 1, 11; 1, 10 lead to a big rip (BR). 
KqF is bounded from below. The numerical solutions marked by (x) are not 
compatible with the main equation l|2.34j) . see Sees. 12.51 and 12.61 for further details. 
Nomenclature: Scale factor a, Hubble scale H = -, acceleration q = event 

* a ' ^ a ' ^^^^ 

horizon radius R and its initial value Rq. For the function KqF see Eqs. (|2.34|) 
and (|2.19jl . for the mass parameters q\,q% see Eqs. (|2.4|) and (|2.6jl . 
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Figure 2.13: The cosmological evolution for different values of today's horizon radius Rq. Here, 

we choose q\ = +2 and q<i = +0, 1. The solutions for Rq = 1, 09; 1, 10 exhibit a big 

crunch (BC), where KqF is unbounded from below. For the initial conditions Rq = 

1,11; 1, 10 the function F and the Hubble scale H approach a finite value, where 

the horizon radius R diverges. The numerical solutions marked by (x) are not 

compatible with the main equation (|2.34j) . see Sees . 12 . 5l and 12 . 6l for further details. 

Nomenclature: Scale factor a, Hubble scale H = -, acceleration g = event 

horizon radius R and its initial value Rq. For the function KqF see Eqs. ([2 .34)1 

and (|2.19[) . for the mass parameters qi,q2 see Eqs. (|2.4j) and ()2.6|) . 
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Table 2.1: This table shows all types of cosmological final stales that we have fuund in Sec. 12.51 
for the all combinations of scaling laws for the CC A (and G for the third law) 
from Sec. 12.21 and the choices for the renormalisation scale fi of Sec. 12.31 Here, dS 
denotes de Sitter solutions, P power-law solutions, BR big rip and BC big crunch 
future singularities, respectively. With exp(±t 2 ) we mean the late-time behaviour of 
the scale factor, where a cx exp(— i 2 ) finally leads to a big crunch. 



2.7 Summary 

Before we summarise the cosmological final states that we have found in the previous 
sections, we have to mention that not all of the solutions are likely to be realised even 
if the underlying assumptions were correct. This is especially important for the extreme 
solutions that exhibit future singularities. For these cases the strength of the gravitational 
field becomes so large that other effects, e.g., higher orders in the curvature scalar or 
unknown quantum gravity effects, cannot be neglected anymore. Also matter sources 
like radiation would become again dominant over dust matter in big crunch scenarios. 
Therefore, it is possible that these cosmological final states are replaced by other solutions. 
For instance, the avoidance of a big rip in phantom cosmologies by such effects has been 
discussed in Refs. jH]- On the other hand, in Ref. [IHj one has an example, where quantum 
effects were not able to prevent a big rip. 

Nevertheless, we stay in this analysis on the basic level and use the results we found 
as an indicator for the (in)stability of the cosmological fate. Indeed, we have found 
regular solutions in many cases, which can be considered as realisable in nature. In this 
sense, we have found de Sitter solutions in the cases of Sees. 12.5. 1| I2.5.2| 12.5.41 and 12.5.81 
Also accelerating and decelerating power-law solutions for the scale factor can be called 
regular, we have found them for all cases with the scaling law ()2.8|) in Sees. 12.5.31 12.5.61 
and l2.5.51 Moreover, we have found in Sec. 12.5.71 a super-exponential expansion law of the 
type a(t) ~ exp(ct 2 ), which however implies negative values of the matter energy density. 
Future singularities of the big rip type are found for all cases with the event horizon as 
renormalisation scale, see Sees. I2.5.4| 12.5.51 and 12.5.61 Big crunch solutions, on the other 
hand, occur only for the scaling law (|2.3jl as described in Sees. 12.5.11 12.5.41 and 12.5.71 In 
Sec. 12.61 we have finally analysed numerically the effect of a running Newton's constant G, 
but we did not observe any changes in the cosmological fates in comparison to the cases 
with constant G. To finish this chapter we give a comprehensive overview of our results 
in Tab. O 



3 Vacuum Energy in Extra Dimensions 



3.1 Introduction 

With this chapter we start to investigate the occurrence of dark energy in models, where 
space-time has some non-trivial properties. One of the most famous approaches for space- 
time modifications represents the introduction of extra spatial dimensions. With the 
concept of Kaluza-Klein (KK) compactification [SO] one has found a framework to unite 
gravity with electrodynamics. Decades later EDs are still a popular way to obtain 4D 
theories from a simpler higher- dimensional setup [51. In this approach, the 4D theory 
which emerges after dimensional reduction is generally characterised by a tower of KK 
modes [221 ■ Also in string theories EDs are a crucial ingredient [S3] and are needed 
for consistency reasons. Therefore EDs have become a common component for modern 
model-building. 

Another interesting aspect of compactified EDs is, that quantum fields in such a non- 
trivial space-time give rise to the Casimir effect [21]. In this scenario the bulk fields have 
to obey certain boundary conditions thereby inducing a finite Casimir energy density, 
which depends on the size and the topology of the EDs. The associated Casimir force 
can be attractive and contract the compactified EDs to a size which is sufficiently small 
so as to have escaped experimental detection so far [231 EH EH] ■ Upon integrating out the 
EDs, the Casimir energies have additionally the interesting property that they appear as 
an effective CC A or vacuum energy in the 4D subspace. In the following chapters we 
will focus mainly on this last point since the Casimir contributions to the vacuum energy 
budget must be smaller than the observed value of the CC in order to avoid fine-tuning. 
One can imagine that this requirement leads to stringent bounds on the size of the ED. 
Later on in Chaps. 0]and[5]we will investigate the Casimir effect for discretised EDs and its 
applications in the context of deconstruction. To prepare for this task we will first discuss 
the calculation the Casimir energy density in the continuum. After that the results can 
be transferred easily to the discretised case. 

3.2 The Casimir Effect 

The Casimir effect is a notable exception from the normal ordering procedure in quantum 
field theories. It occurs when quantum fields have to obey certain boundary conditions, 
for instance the electric component of the photon field, restricted between two parallel 
conducting plates, has to vanish on the plates. This causes a geometry dependent vacuum 
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energy density inducing a force on the plates. Therefore, the Casimir effect is a macro- 
scopic quantum phenomenon, which is experimentally well established [SB]- For some 
recent reviews of the effect and its applications see Refs. (211 EE]- Global properties of a 
non-trivial space-time may also be probed by the Casimir effect which is sensitive to the 
IR structure of a theory. 

Let us now consider the Casimir effect for a scalar field in a 5D theory, where the ED 
is compactified on the circle S 1 with circumference R. Hence, the whole manifold has the 
topology M. x S 1 , where Ai denotes Minkowski space. In order to find the possible field 
theory configurations we have to specify the boundary conditions for the scalar field on 
the circle. Motivated by periodicity one obvious choice is given by 

4>{y + R) = 4>{y), (3-1) 

where y is the coordinate in the fifth dimension. This choice implies a cylinder-like 
structure. But this not the only one. According to Ref. [HE] one can also form a structure 
with the properties of a Mobius band, where <fi obeys anti-periodic boundary conditions 

<f>(y + mR) = (-l) m <f>(y), m G Z, (3.2) 

In this case one must cycle twice through the circle S l to completely traverse the Mobius 
band. In the latter case, the field is called a twisted field, whereas the fields with periodic 
boundary conditions are called untwisted fields. Locally, both cases have the same product 
structure, but globally they differ significantly. Since they yield inequivalent degrees of 
freedom of the field 0, both must be considered in the Casimir effect. 

Turning to the calculation of the Casimir energy, we call q the momentum corresponding 
to the position y in the extra-dimensional space. Since the 5D manifold is flat, it is sensible 
to use for the scalar field a plane wave ansatz given by 

4>(t, x,y) — J ■ exp(iut — ipx — iqy), (3.3) 

where J is a normalisation factor, uo the energy and p the 3-momentum corresponding to 
the time t and respectively the 3-coordinates x. As discussed above, the untwisted field 
configuration is fixed by periodic boundary conditions, 

4>i3l + R)=m =► e"^ = l, 
implying a discrete momentum spectrum, 

71 

q = 2tt-, n£Z. (3.4) 
R 

For twisted fields we have to use anti-periodic boundary conditions 

(f)(y + mR) = (-l) m 0(y), mGZ e~ iqR = -1, 



which yield the discrete momentum spectrum 
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Since this is the only difference which is relevant for the following calculations, we will 
work with untwisted fields and replace n by n — 1/2 when needed. 

To normalise the field modes in Eq. ()3.3)1 . we define the following scalar product for two 
modes </>i i2 (p, n) by 

(01,^2) := i Jd 3 x j dy [f x {d t (j> 2 ) - ((^OV) , 

so that the normalisation factor J can be fixed by demanding the orthonormality relation 

(0(p,n),0(p>')) = -Vi l 5®(p-p')5 nn ,, (3.6) 

where V3 is an arbitrary 3-volume factor, which leaves the scalar product dimensionless. 
With the ansatz in Eq. ()3.3j) . we find 

jtj= MW' (3 ' T) 

where we have applied the relations 

/f R n — n' 

d 3 x • exp(ix(p — p')) = (2ir) 3 S ( ' 3 \p — p') and / dy exp(2iriy — - — ) = R8 nn /. 
Jo R 



The equation of motion for the real 5D scalar field with a bulk mass M s is given by the 
Klein-Gordon equation 

dt 2 



- V 2 - 4d 5 + M s 2 



= 0, (3.8) 
which determines the energy u of a field mode with the momenta p and q: 



2 ->2 , 2 

u = p + m . 
Here we have introduced the squared effective 4D mass 

m 2 := q 2 + M s 2 . 

Furthermore, we need the 5D energy-momentum tensor Tab of the real scalar field <fi to 
calculate the energy density and pressure of the field. It has the form [TS] 

Tab = (0,A) f (0,B) - \gAB9 CD {<P,c)K<P, D ) + l -g AB M 2 s ^^ (3.9) 

where A, B,C, . . . are 5D coordinate indices. Here, A = is the time-like index, A = 1,2,3 
are spatial indices, corresponding to the uncompactified 3-space, and A = 5 characterises 
the extra spatial dimension. The 5D energy density ps is then given by the 00-component 
of Tab- 

p 5 = Too = hdo^iM) + i(V0) f (V0) + ^(9 5 0) t (9 5 0) + ^M 2 0t0. (3.10) 
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By averaging over all directions of the isotropic 3-space, we obtain the pressure p^ of the 
scalar field 6: 



1^ 1 



8=1 



(V0)t(V0) + ^(<9 o 0)t(<9 o 0) - \^<t>)K^<t>) - ^(^W) - 

(3.11)' 

Let us now perform the canonical quantisation of the field <fi by introducing the field 
operator 

r N 

(3.12) 



y) = Vvs / d 3 p (^(p> ra ) fl p> + + (p> ™)4>) ' 

n=l 

where a Pyn and a£ n obey bosonic commutator relations: 

*P,n> a l>,n>) = V Z 1 5{P- P')$nn> , [«, «] = [a 1 , « f ] = 0. 



Now, the 00-component T 00 and the averaged ii-components ~ XliLi ^ °f the energy- 



momentum operator Tab follow from substituting the field operator in Eq. (|3.12|) into 
Eqs. ()3.10p and ()3.11|) . Here, it is useful to consider the relations 

p p oo 

{d^)\d^) = Vl J d 3 p J d 3 p' [ d ^& n K> ■ <W, n')% tn , + ■■■' 

n,n'=—oo 

/oo 
d 3 p ^A-u 2 + --- 
n=—oo 

/oo 
d 3 p A ] A-p 2 b + --- V6= 1,2,3 

n=— oo 
» oo 

V 3 J d 3 p ^A-q 2 + ---, (3.13) 



where the ellipses (•■•) denote the terms which vanish in the VEVs (0|Tab|0) due 
to (0\aJ = a\0) = 0. When we insert the terms from Eqs. ()3.13j) into Eqs. ({HUD} 
and ()3.11|) . one obtains, after taking the VEVs of T 00 and T#, the energy density p 5 
and the pressure p 5 of the quantised field 0, 



P5 



/oo 
d 3 P J]J - 
n=—oo 

/oo 
d 3 p J]J - 
n=—oo 

1 3 „ oo 

3E<°i f -i°> = vsl d 3 P J2 

i=l " n=—oo 



-p 2 + -uo 2 - -p 2 - -q 2 - -M 2 
3^ 2 T 2 q 2 s 



1 ■/••'>!: 



3d; ' 
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where we have used the energy-momentum relation uj 2 = p 2 + q 2 + M 2 . With the normal- 
isation factor J from Eq. (|3.7|) we finally obtain the energy density p 5 and the pressure p 5 
of the quantised 5D field 0, 



1 f °° 

/^=(0|:/nn|0) = -7—777 / d"7' £ UJ, (3.14) 



P5 

i=l 



2(2tt) ; 



n=— oo 

3 ., , oo 



The momentum integral j d 3 p = Ait J °° dp • p 2 in these equations and thus ps and ps are 
divergent. We therefore have to apply a regularisation procedure to obtain meaningful, 
finite expressions. Consider first Eq. f)3. 14jl . By introducing an exponential suppression 
factor e~( dp ^ as regulator function for d > 0, we find 



dp • p 2 ui — ► / dp • p 2 cue 



-(dp) 2 

= ^rrtd^eS^K^ (\d 2 m 2 ) (3.16) 

= -cT 4 + -m 2 d~ 2 + -m 4 ( - + -7 - In 2 + \n{dm) j + C(rf 6 m 6 ), 
2 4 8 \4 2 / 

where u = \Jp 2 + m 2 and K n (x) denotes the modified Bessel function of the second kind 
of the order n and 7 = 0,577. . . is the Euler-Mascheroni constant. The limit d — > 
removes the regulator and recovers the divergence. Before taking this limit, a renormali- 
sation has to be carried out to remove the potentially divergent terms. Alternatively to the 
exponential regulator function in Eq. ()3.16|) . one can also apply dimensional regularisation 
by moving to n space-time dimensions. To be specific, 

dp • p (n " 2) • UJ 
2 \mJ V 2 



n— >4 / . x 1 4 



(47r)-m 4 



(;/ ' • I" ( ^ ) +0(n - I) 



1 

- n^4 



where the renormalisation scale p has been introduced to keep the mass dimension of the 
whole term constant, and S(n) is the surface area of an n-ball. The regularisation of the 
divergent integral in Eq. ()3.15j) with the exponential suppression factor with d > goes 
along the same lines as above: 

Jl POO 9 



1 2P 

dp • p — 

UJ 



poo 2 

/ dp .p2^ e -(*) 2 (317) 
Jo w 

l -m 2 d- 2 e^ 2 (d 2 m 2 K {\d 2 m 2 ) + (1 - d 2 m 2 )K 1 {\d 2 m 2 )) 
^cT 4 - X -m 2 d- 2 - jjm 4 (J^ + \ ~ In 2 + \n(dm)j + 0(d 6 m 6 ). 
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In the next step we have to remove the potential divergences in the regularised expres- 
sions. For a curved background space-time [TH], one would therefore decompose p into a 
divergent and a finite term, so that the former one has the form of a cosmological term 
in Einsteins's equations. Then the divergences would be absorbed to yield renormalised 
coupling constants (which is the origin of some RGEs for A and G in Chap. EJ, and the 
finite remainder is called the renormalised energy density or, in our case, the Casimir 
energy density. Here, such a general treatment is not necessary because the divergence 
also arises in flat space-time, like our Ai x 5 1 -manifold, but there are neither cosmolog- 
ical terms nor Einstein's equations. In order to get rid of the divergence, one simply 
subtracts the corresponding part of the energy density of the same field in a Minkowski- 
like space-time with the same dimensions, i.e., M. x IR 1 in our case. In this case there 
is no IR cutoff and one has to integrate over the 5-momentum instead of summing over 
it. This kind of renormalisation works because the 5D Minkowski space suffers from the 
same divergence as the A4 x S 1 space-time but exhibits no Casimir effect. For the energy 
density ps from Eq. (|3.14jl . where the p-integral has been regularised by using Eq. (j3.16|) . 
the renormalisation subtraction can be written as 



P5,renorm OC 



00 POO 

J2 f(ji) - / dn/(n) (3.18) 
„— ~, J —00 



with 



f{n) ■= IcT 4 + im 2 cT 2 + ^m 4 Q + ^7 - In 2 + ln(rfm)^ + 0(d 6 m 6 ). 
Above we have applied the substitution q = 2ixn/R to rewrite the integral over q: 



2tt j y1 ' RJ J y R 

Notice that when keeping in Eqs. ()3.16|) and ()3.17|) only terms proportional to m 4 lnm, 
we obtain the equation of state p = —p of the CC. As we will see next, all other terms 
including the divergences will vanish due to the subtraction in Eq. ()3.18)1 and the sub- 
sequent regularisation removal with d — > 0. The finite result of the subtraction can be 
calculated by using the Abel-Plana formulas jHO] given by 



E/(») 



n=0 



n=0 



„ + !)_, in . Kn) = ->l d „. Zm(zM . ( 3. 20 ) 



2 

m , 



With m = ^Jk 2 n 2 + M s 2 and k := 2tc/R we immediately see that terms like mF 
m 4 in the function f(n) are cancelled on the right-hand side. And the term |/(0) can 
be dropped as well since we consider the range n = —00 ... 00. Only terms of the form 
m 4 In m survive, which we write for untwisted fields in the form 

A := m 4 (n) ln[m(n)] — / dn • m 4 (n) ln[m(n)]. 
, J~ 00 
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By applying the first Abel-Plana formula (J3.19)) we therefore obtain 

m 4 (+m) ln[m(+m)] — m 4 (— in) ln[m(— in)] 



A = 2i dn 



o 



exp(27m) — 1 



Assuming k,n,M > 0, we have to consider two cases for the root m(±in) in the last 
equation: 

(M 2 -k 2 n 2 )^ for M > kn, 
±i(k 2 n 2 - M 2 )i for M < kn. 



[k 2 (±tn) 2 + M 2 \ 



For a > we can write the logarithm as ln(±i • a) = ±m/2 + In a, and with x := M/k we 
obtain the result 

(n 2 -x 2 ) 2 



A(M) = -2nk* I dn^- ^ — , (3.21) 

' exp(27rn)-l' 



which has, in the massless case (x = 0), the value 

A = -27^^- C(5) = -8iT 4 3C(5). 

For twisted fields we have to use the second Abel-Plana formula ()3.20|) . Analogously, we 
write 

f°° m 4 (+m) ln[m(+w)] — m 4 (— in) ln[m(— in)] 

B(M) = —2i I dn — — r , 

V ; J exp(2vm) + l 



using m(n) = \/k 2 n 2 + M 2 , and not m(n) = y k 2 (n + \) 2 + M 2 . Thus, the result 
becomes 

B(M) = +2nk 4 / dn^- (3.22) 

v ; J x exp(27rn) + l' v ; 

which has for massless fields (x = 0) the value 

B = +27rA;4 64^- C(5) = Y i? " 43C(5) - 

For large masses (x ^> 1), approximate expressions for A(M) and B(M) can be given by 
neglecting the 1 in the denominator of the Eqs. (j3.21|) and ()3.22|) : 

exp(27rn) ±1 exp(27rn) 47r 5 

Combining the prefactors in Eq. ()3.14j) and the finite results A(M), B(M) from the renor- 
malisation in Eq. (j3.18|) we have found the finite 5D energy density. To obtain the effective 
4D Casimir energy p 4 density we just have to integrate over the fifth dimension j Q R dy, 
which simply yields a factor R. Finally, we find for untwisted scalar fields 

P4 = ^y 2 - A (M) M =° -1 • (2vrr 2 /T 4 3C(5), (3.24) 
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whereas twisted scalar fields yield 

P4 = ■ B(M) M =° +^ ■ (2vr)- 2 i?r 4 3C(5). (3.25) 

Obviously, untwisted and twisted fields provide energy densities p 4 of different sign. Note 
that the values for p 4 in Eqs. (13.24(1 and (|3.25(1 agree with the results in Refs. [^3 Ell 
for M = 0. 

As explained above, the corresponding effective 4D pressure p 4 is just the negative of p 4 . 
Therefore the Casimir energy density from the ED represents a finite contribution to the 
4D CC. From these results we see that p 4 scales like the inverse fourth power of the size R 
of the fifth dimension. For small EDs this could cause problems with the observed tiny 
value of the CC. Apart from fine-tuning the field content, there is another possibility 
to solve this problem. From expression (|3.23(1 we obtain an approximately exponential 
suppression of the Casimir energy by bulk field masses M = xk. We will make use of this 
behaviour in Sees. 14.51 and 15.21 where the Casimir energy becomes small for large bulk 
masses even when the ED is small. 

Finally, we mention that from statistical arguments and counting degrees of freedom 
one can immediately conclude that the Casimir energy densities for Dirac fermions are 
just (—4) times the values of real scalars. We will use this later in Sec. 14.31 



4 Discretised Extra Dimensions 



4.1 Introduction 

In the Chap.Elwe have discussed dark energy in the form of Casimir energy in the context 
of a continuous ED. We have found there that the Casimir effect yields a contribution 
the effective 4D CC. In this chapter we go one step further and consider discretised EDs. 
As already mentioned in Chap. ^ a discrete space-time structure might not only serve 
as an UV regulator in theories of quantum gravity but it is also a crucial ingredient of 
deconstruction. In this framework the phenomenology of discretised higher dimensions 
is exactly reproduced by a deconstruction model working in a 4D continuous space-time. 
This concept helps to circumvent some disadvantages occurring in extra-dimensional the- 
ories. 

In the following sections we will derive the Casimir energy density for scalar and fermion 
fields in a discretised ED. Since we we are working effectively on a transverse lattice 
we expect to encounter typical lattice effects like fermion doubling. Furthermore, the 
scaling of the Casmir energy density with the number of lattice sites and the suppression 
by bulk masses will be discussed in detail. In Chap. El we will apply these results to a 
deconstruction model and propose a prescription to calculate the absolute value of the 
vacuum energy of 4D quantum fields in the deconstruction framework. 



4.2 Casimir Effect for a Scalar Field 

In this section, we consider a scalar quantum field in a space-time with the topol- 
ogy M. x «Sj^ t , where M. is the continuous Minkowski space and <S^ t denotes the discrete 
fifth dimension compactified on the circle. Taking the discrete nature of the fifth dimen- 
sion into account, the discretisation of the circle S 1 also forces the coordinate y in the 
fifth dimension to be discrete. Assuming N lattice sites with a universal lattice spacing a, 
the circumference of the fifth dimension is given by R = Na, and the position y of each 
site can be described by a coordinate index j, 

y = a-j, j = l,...,N. (4.1) 

From the standard definition for a derivative in the continuum, 

^(y) = lim^ + a) "^ ) > 
oy a-*o a 
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follows the discrete forward and backward difference operators d§(f) and (<9s0)^: 

ay a 

By inserting the ansatz (|3.3J) for we find 

<9 5 = oT\e-^ a - 1)0, (950) 1 " = a~\\ - e^)0, 

and therefore 

<95<9 5 = — 2a~ 2 (l — cosga)0. 

Taking into account the discrete derivatives and c^, the continuum Klein-Gordon equa- 
tion (|3.8J) for a real 5D scalar field with bulk mass M s becomes 



0. 



Thus the the energy uj of a field mode with the momenta p and g reads 

u 2 =p 2 + m 2 , m 2 := 2cT 2 (l - cosga) + M 2 , (4.2) 

where the values g depend on the boundary conditions as in Eqs. ()3.4j) and (|3.5|) . However, 
here n has to be in the range 1 ... AT since the lattice introduces an UV cutoff. Values 
of n outside the range 1 . . . N are mapped back into this range due to the 27r-periodicity 
in the discrete derivatives. Therefore, the momentum spectrum is finite and given by 

g = 2 7r£ = ^.£, n = l...N (4.3) 
R a IS 

for untwisted fields. For twisted fields we respectively find 

(n - |) 27T (n - 1) , . 

q = 2n K -—^ = — K -—^, n = l...N. 4.4 

R a IN 

We can now easily transfer many of the continuum results from Sec. 13.21 to the discretised 
case by replacing the infinite sum Y^=-oo by the finhe sum Yun=v For instance, the 
scalar product for the field modes now reads 

f N 

(0 1; 2 ) := i d 3 xaJ2 U\(fr<h) ~ (^0i) f 02 
J i=i 

which yields via Eq. ()3.6j) the normalisation factor J from Eq. (|3.7J) by using the relation 

N 



^2 ex P ( 2ni —jf-j ) = N5 nn'- 
3=1 ^ " ' 
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The quantisation of the scalar field and the determination of the corresponding energy 
density p 5 and pressure p 5 goes along the same lines as in the continuum case: 



P5 = (0|T 00 |0) 



1 



2(27r) 3 i? 

P6 = £ X>| f« |0) ' 



N 



n=l 



i=l 



2(2n) 3 R 



N _n 

1 ^ 3u 



(4.5) 
(4.6) 



n=l 



Also, the regularisation of the momentum integral J d 3 p can be performed exactly like in 
the continuum case. Here, we also consider the method with the exponential regulator 
function e~^ p ^ with d > 0. However, in the renormalisation procedure we cannot make 
use of the Abel-Plana formulas since they only treat infinite sums. In order to obtain the 
finite Casimir energy density, it is necessary to compare the discrete mode sums belonging 
to the momenta in S^ t with the energy density and pressure of a field in a space-time with 
a non-compactified but still discretised ED. Regarding Eqs. ()4.5j) and (|4.6jl . the mode sum 
with respect to the fifth momentum coordinate q = 27m/(aN) is of the type 



N 



£/(n/A0, 



n=l 



where f(n/N) summarises the terms 



+ l m 2 d~ 2 + -m 4 



i + l 7 -ln2 + ln(ci)J 



0(d 6 m^ 



following from the last line in Eq. ([3. 16)1 . From this sum, the mode integral corresponding 
to a non-compactified ^-dimension can be obtained by cutting out a section of length R 
of an R 1 -dimension. This means, that we take the limit of an infinite number M of lattice 
sites, M — > oo, while keeping the spacing a constant: 



R_ 



M 



n=l 



n 
M 



R 



' M A 

y—f (■ 

^ M J \M 

71=1 



n 



=n/M 



M^oo 



N- / ds-f(s), 



where Ma becomes the infinite "length" of M 1 and f(s) is the same function as in the <S^ t 
mode sum. In the last equation, we have substituted s := n/M and inserted An = 1 
so that ds = An/M for M — > oo. Both the sum and the integral are finite since the 
lattice introduces an UV cutoff. Then the renormalisation is performed by subtracting 
the integral from the sum, 



N „i N i 

J2f(—)~ N - ds- f(s) = ym 4 \nm- N ■ ds-m 4 ln 

71=1 ^° 71=1 ^° 



m, 



(4.7) 



where only m 4 lnm survives since all other terms like m°, m 2 , m 4 either vanish when the 
regularisation is removed for d —>■ or are completely subtracted due to the following 
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identities: 



N 

D 

n=l 



l-cos2vrf + |a 2 M s 2 ) 



JV ■ I ds ■ (1 - cos 2tts + ±a 2 Mf 
'o 



= jv(i + Wm*), 



(4.8) 
(4.9) 



V(l-cos27rf + ±a 2 M 2 ) 2 = JV • / ds ■ (1 - cos 2vrs + ±a 2 M 2 ) 2 (4.10) 

n=l " 

= JV(f + a 2 M 2 + (±a 2 M 2 ) 2 ). (4.11) 

This is also the case for twisted fields, where n is replaced by n — \. Note that in the 
discretised case we obtain the CC equation of state p = — p, too. Finally, we find the 
renormalised Casimir energy density to be 



1 



P5 



47T 



N ! 

> m 4 lnm- JV • / ds-m 4 lnm 



2(271)3/? 8 

= +iT 5 1 Si(JV) = -p 5 , (untwisted) (4.12) 
where m 2 = 2a _2 (l — cosga) + M 2 and where we have introduced the function 

N / N 2 



Si(iV) 



1 



4(2tt) 



r JV 4 - 



E (i - cos27r ^ + ^ 2M s 2 ) ^ (l - cos2tt^ + \a 2 M^j 

. (4.13) 



A' ■ / ds- ( 1 -cos2tts + -a 2 M s 2 j In ( 1 — cos 2ns + -a 2 M 2 



In the limit JV — > oo and M s = 0, the function Si(N) converges to the value of a continuous 
fifth dimension: 

lim 5l (JV) = --i--3C(5).4. 
By integrating out the fifth dimension, we obtain the 4D energy density 

Pa = £ dr • p 5 = R ■ p 5 = ~ = Jf4 * ("°> 0787970 • • • )■ (untwisted) (4.14) 

In the case of a twisted scalar field everything is like above, but the energy density reads 

P5 = +ir 5 s 2 (jv) 

= — p 5 , (twisted) 

where S 2 (N) is the function Si(JV) with n replaced by n — \ . For massless fields (M s = 0) 
we obtain in the continuum limit 



] im .^(iV) = + 77 ^3C(5)-(4-i), 



N^oo 



4(2tt) 2 
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and after integrating out the fifth dimension the 4D energy density reads 

15 3C(5) 1 

p 4 = R.p 5 = +-. — ^ = — ■ (+0, 0738722 . . . ). (twisted) (4.15) 

In the continuum limit the lattice results for p 4 found here are consistent with the con- 
tinuum Casimir energy densities from Eqs. 1)3.24)1 and (|3.25j) . 



4.3 Casimir Effect for a Dirac Fermion 

In analogy with the treatment of scalar fields in previous section, we will now calculate 
the Casimir energy density of Dirac fermions. Therefore, a plane wave Ansatz for Dirac 
spinor fields \l/ in the M. x S 1 manifold is a convenient choice, too: 

^ = ?/>exp(— iujt + ipx + iqy). (4-16) 

The boundary conditions, associated with the compactified ^-dimension, provide the 
discrete momentum spectra. For twisted and untwisted fields we have as before q = 
2n(n — |)/(aiV) and q = 2nn/(aN), respectively. Like in Eq. ()4.1)) the coordinate y 
corresponding to the fifth dimension is discrete, y = a- j, where j = 1, . . . , N, and implies 
an upper bound for the momentum q. 

Unlike the Klein-Gordon equation for scalars fields, the Dirac equation is linear in the 
derivatives, and therefore we need a symmetric derivative operator for the discrete y- 
coordinate: 

With the Ansatz ()4.16j) we obtain 

d§fy(j) = —ip ■ exp(— iojt + ipx) [exp(iqa(j + 1)) — exp(iqa(j — 1))] 
= ■ (V-Uin(ga)^ , 

and together with the 5D Dirac equation 1 for a Dirac field with mass Mf, 

(i>y A d A -M { )* = 0, A = 0,...,3,5, 
the energy-momentum relation is determined to be 

u 2 =p 2 + m 2 , m 2 := a~ 2 sin 2 qa + Mf . (4.17) 
The energy-momentum tensor Tab f° r the Dirac field $ has the form JE] 

Tab = -i[*7A^ B * + *7b<9a* - («9 a ^) 7b ^ - {d B V)i A V], 

1 A fifth Dirac matrix 7 s := ijl D has to be introduced, where 7| D is the usual 7 s matrix of the 4D Dirac 
theory 
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and the usual canonical quantisation procedure parallels that for scalar fields up to replac- 
ing the bosonic commutator relations by the fermionic anti-commutator relations, which 
give an overall minus sign in the result. The Dirac fermion also has four times the degrees 
of freedoms of a real scalar, describing particles and anti-particles with two spin states 
each. In total, the energy density ps and pressure ps of a quantised Dirac field differ from 
the scalar results of Eqs. ()4.5)1 and ()4.6)1 only by a factor of (—4) and in the modified 
energy- momentum relation of Eq. f)4 . 1 7|) . i.e., 



P5 



N 



v ' 71=1 



p-u, 



o N r ->2 

2 E/dV 



CO 



(4.18) 
(4.19) 



where u 2 

Eoo 
n- 



p 2 + a 2 sin 2 qa + M 2 



By replacing back the sum J2n=i ^° ^ ne continuum 



and additionally using the continuum energy momentum relation 



p 2 + q 2 + M 2 , we now explicitly see by comparison with Eqs. ()3.14j) and (|3.15jl that the 
corresponding fermionic Casimir energy density would indeed be (—4) times the scalar 
values. In the following, however, we will observe that on the lattice this is not true 
anymore due to fermion doublers. 

From here on, the regularisation and renormalisation procedures are identical to the scalar 
case in Sec. 14.21 This also implies that the equation of state of the fermionic vacuum energy 
is that of a cosmo logical constant, p = —p. Thus, it is sufficient to give the renormalised 
energy density in five dimensions 



P5 



4:71- 



{2TifR 
+R~ 5 F 1 (N), 



N i 

>^ m 4 In m — N ■ / ds ■ m 4 In 
tl Jo 



m 



(4.20) 



where the function F\(N) in the last equation is defined as 

N 



Fi(N) 



1 



4(2vr) : 



-N 



( sin2 27T jj + a 2 M? ) In (sin 2 2vr^ + a 2 M f 2 

n=l 



N ■ ds - ( sin 2 2vr— + a 2 M f 2 ) In I sin 2 2tt— + a 2 M f 2 



with R = Na. For the twisted Dirac field we have 



(4.21) 



p 5 = R~ 5 F 2 (N), 

where F2(N) is the function Fi(N) with n replaced by n — \. Unlike the functions Si^(N) 
for the scalar fields, the functions F 12 (iV) for the fermionic fields have two limit points 
each, which depend on whether the number of lattice sites N is even or odd. For massless 
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fermions (Mf = 0) and even iV we obtain 

1 



lim Fx (AH 
lim F 2 (N) 



■3C(5) ■ (-32), (untwisted) 



4(2tt) 2 

-— ^ ■ 3C(5) • (+30). (twisted) 



After integrating out the fifth dimension, the 4D Casimir energy densities read 



Pi 



Pi 



32 • 3((5) 
4(2vr) 2 J R 4 " 
-30 ■ 3C(5) 



— •(+0,630376 



1 

B 1 



;-0, 590978. 



4(2tt) 2 # 4 

In the case of odd A", both functions have the same limit 

lim FAN) = lim F 2 (N) = — -5— 
w-oo v 7 jv-^oo v 7 4(2vr) : 



(untwisted) 
. (twisted) 



(4.22) 
(4.23) 



3C(5), 



Pi 



3C(5) --L- (+0,019699... 
it 



4(27r) 2 i? 4 

Obviously, this behaviour is an effect of the lattice, in Ref. it is called an odd-even 
artefact. We also notice, that the N —> oo limit of the sum of twisted and untwisted results 
does not depend on whether A" is even or odd. Therefore, it is reasonable to consider only 
this sum as a physical quantity. For finite N, this odd-even artefact is illustrated in 
Fig. 14.11 Note again, that the continuum results for p 4 in Eqs. (|4.22j) and ()4.23|) are 
identical with the values in Ref. 



4.4 Massless Fields 

The calculations of Sec. 14.21 show that the Casimir effect for a real scalar field in the 
transverse lattice space-time M. x S\ &t induces a negative vacuum energy density p and 
therefore a negative contribution to the effective 4D CC. On the other hand, the fermionic 
Dirac field of Sec. 14.31 yields a positive contribution to the CC. 

We have already concluded that only the sum of twisted and untwisted fields can be 
regarded as a physical quantity, and we note that its sign is independent of N. More- 
over, for a constant circumference R and small N, the Casimir energy density Pa{N) in 
the transverse lattice setup has already the same order of magnitude as the energy den- 
sity Pa(N — ► oo) in the continuum limit. Specifically, for A^ > 10 the continuum result is 
approximated at the few percent level. Even for a number of lattice sites which is as small 
as A" = 3, the results differ at most by a factor of 2, which is clearly shown in Fig. 14.21 

In the limit N oo (Tab. 14. 1|) . the results for real scalars are the same as in the non- 
lattice calculation [23 but for the fermions there is an extra factor of 2 in the energy 
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Figure 4.1: Illustration of the odd-even artefact for fermion fields. The cases of odd and even N 
are plotted separately, and the circumference R of the fifth dimension is kept con- 
stant. For a better representation, we flipped the sign in the results for twisted 
fields and even N. The horizontal lines correspond to the continuum values given 
in Table O 



density of our lattice calculation because of the fermion doubling phenomenon in lattice 
theory. In a calculation for continuous dimensions, one usually expects, from counting 
degrees of freedom, that the energy density for Dirac fermions is (—4) times the value of 
real scalars. 



Up to now, we have investigated the Casimir effect for Dirac fermions and real scalars 
having twisted and untwisted field configurations. When passing to a complex scalar 
field which transforms under a U(l) gauge group there exist only trivial (untwisted) 
structures and therefore the charged scalar obeys only periodic boundary conditions. For 
fermions, on the other hand, the appearance of twisted field modes is related to the double 
covering map SL(2, C) — > 50(3, 1) which gives rise to inequivalent spin connections [59] . 
Consequently, even in presence of a simply connected gauge group like U(l) we still have 
also the anti-periodic boundary condition for the fermions. 
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Figure 4.2: The Casimir energy densities for fermions (upper graph) and real scalars (lower 
graph). In this plot, R is fixed so that the lattice spacing decreases for increasing N. 
The dashed horizontal lines denote the continuum limit N — > 00, which has for 
scalars the value — jg and for fermions +5. 
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^'exponential suppression factors iti Kqs. (|4.24j) ana (|4.25|) ot the Casrrrnr energy 
densities for TV = 3 lattice sites, where M Sj f denotes the bulk masses of the quantum 
fields. 

4.5 Exponential Suppression by Massive Fields 



So far, we have given results only in the case of vanishing bulk mass (M B) f = 0). For 
massive 5D fields we observe an approximately exponential suppression of the Casimir 
energy This behaviour becomes obvious in the analytical calculation for a continuous 
ED, which is given in Sec. 13.21 But it is also achieved for the discretised case of this 
chapter, where in the limit of an infinite number of lattice sites (N — > oo) we approach 
the values of the analytical formulas (|3.24J) and ()3.25j) . To investigate the suppression 
behaviour depending on the mass M S) f and the number N, we examine the ratio between 
the energy density of fields with mass M s f and that of massless fields. For scalar fields 
this ratio is defined by 

S^MM) + S 2 (M S R) 

s 1 (o) + s 2 (o) ' 1 • ] 

where the functions Si^{N) are taken from Eq. ()4.13|) of Sec. 14.21 Analogously, using the 
functions Fi 2 (N) from Eq. (j4.21|) of Sec. 14.31 the ratio for fermionic fields reads 



FjjMtR) + F 2 (M{R) 
F 1 (0)+F 2 (0) 



(4.25) 



Both ratios are plotted in Fig. 14.31 for a range of values of iV and M s jR. The suppression 
by a bulk mass is most minimal for small N. In the case of iV = 3 lattice sites, the 
corresponding ratios are given in Table 14.21 
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Scalar Fermion 




MR 

Figure 4.3: The ratios in Eq's. (|4.24() and (|4.25f) of the Casimir energy density between massive 
and massless fields. The values for N — > oo are taken from the analytic formulas in 
Eqs. (HQS) and (ET251) . 



5 Vacuum Energy in Deconstruction 



In this chapter we will discuss the vacuum energy of 4D quantum fields occurring in a 
deconstruction scenario. It is well known that the absolute value of the zero-point energy 
density J d 3 p ■ cu(p) of 4D quantum fields is formally divergent and completely uncon- 
strained from a theory point of view. This means that without any further information 
we do not know how to assign a sensible and finite vacuum energy value to these fields. 
In Sec. 12.11 we already mentioned this problem and investigated the change of the vac- 
uum energy density with respect to a renormalisation scale, but the absolute value still 
remained undetermined there. 

Here, however, we propose a prescription how to treat this problem by using the corre- 
spondence of a discretised fifth dimension and deconstruction. For a 5D field obeying 
certain boundary conditions we are able to derive the part of its vacuum energy that 
depends on the boundary via the Casimir effect. Since in deconstruction the 5D field is 
described by many 4D fields it seems plausible to identify the vacuum energies of all the 
4D fields with the Casimir energy density of the 5D field. Going this way we have found a 
definite and sensible prescription to determine the zero-point energies of quantum fields in 
the deconstruction framework. Due to the lack of alternatives to handle this problem, our 
idea represents at least a well motivated ansatz to obtain a reasonable result, in contrast 
to the naive cutoff method given in Eq. (|2.1jl . 

In the next section we will briefly introduce a deconstruction model that describes a 
discretised fifth dimension and give the values of masses and VEVs of the fields that 
appear in this framework. After that, we will use the results of the previous chapter to 
determine the vacuum energy of the 4D fields by using the prescription given above. 

5.1 Deconstruction Model 

In continuous EDs the maximum number of KK modes is usually restricted by an UV 
cutoff which reflects the fact that non-Abelian gauge theories in higher dimensions are 
non-renormalisable. Although this leads below the cutoff to a renormalisable effective 
4D theory, the full higher-dimensional gauge- invariance is in general lost. To circum- 
vent these problems the deconstruction framework has been proposed [23 EDI represent- 
ing manifestly gauge-invariant and renormalisable 4D gauge theories, which reproduce 
higher-dimensional physics in their IR limit. These theories use the transverse lattice 
technique [61] as a gauge-invariant regulator to describe the EDs and yield viable UV 
completions extra-dimensional theories f65| . 
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4> ■ St '4-1 St '4-9 

Figure 5.1: Part of the moose 'diagram, where the sites i, i + 1 and i + 2 are shown for a 
deconstructed extra dimension compactified on S . Each circle Gi corresponds to the 
U(l)i gauge group and each arrow pointing towards (outwards) a circle represents 
a field with negative (positive) charge under this group |66j . 

Here, we consider a specific 5D deconstruction framework which was introduced in Ref. [HHl 
and which represents a periodic model for a deconstructed 5D U(l) gauge theory com- 
pactified on the circle S 1 . The model is denned by an U(1) N = Il^ 1 C/(l)j product gauge 
group with iV scalar link fields Qi (i = 1, . . . , N), which carry the [/(l)-charges (q, —q) 
under the neighbouring groups U(l)i x U{l) i+ i. The identification i + N = i establishes 
the periodicity of the lattice corresponding to untwisted quantum fields in the language 
of Sec. 13 .2\ whereas twisted quantum fields are described by an anti-periodic lattice with 
the condition i + N = —i. On the i th lattice site, we put one Dirac fermion and 
one scalar <3>j which carry both the charge — q of the group Z/ The fermions are 
SM-singlets and correspond to a right-handed bulk neutrino in the ADD scheme jHZl EE] ■ 
As an illustration of this setup the corresponding "moose" |HH] (or "quiver" [7U] ) diagram 
is shown in Fig. 15.11 Let us split up the Lagrangian of the model into several parts, 

£ = £ kin [$ 4 , + C kin [A^] + £ kin [^] + C mass [^ h Qi] - V, 

where £kin[^4f] and £kin[^i] are the standard kinetic terms for the gauge bosons A? and 
the fermions given by 

1 N 

£\dn[Ai] = —— ^J(^Aii/ — d u A nfl ) 2 , 

71=1 

N 

Adn[^] = ^^nl^-ignA^n. 
n=l 

Moreover, the kinetic terms /^m^i, Qi] for the scalars $j and Qi, which will later provide 
the gauge boson masses, can be written as 

N 

|(^u +igiA ili )$ i \ +\{dp + igiAi^-ig i+1 A { i + i)i 1 )Qi\ , (5.1) 

t=i 

where <?j are the gauge couplings corresponding to the gauge group Z7 Finally, the 
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mass and mixing terms involving the fermions and the link fields Qi are combined into 



U 

2 



N 



n=l 



■*(n+l)R *(n-l)R 



+ h.c., 



(5.2) 



where ^l,r := |(1 T 7 5 ) denote the left- and right-handed components of the Dirac 
fermion In the last line we have already introduced the universal real VEV (Qi) = u 
of the link fields Qi. At this point let us denote the VEV of the scalars by = v, 
which is also real and universal for alH = 1 ... N. These properties of the VEV structure 
follow from the renormalisable potential V, which is not shown here since it not needed for 
the determination of vacuum energy in this setup. The potential has been discussed and 
minimised in great detail in Ref. [HE], where it leads via a type-II seesaw mechanism [71] 
to the following values for the VEVs: 



(Qi) = u~ 10"' eV, = v ~ 1(T GeV. 



Once Qi and have acquired their VEVs the gauge bosons A ifl become massive via the 
Higgs mechanism. Assuming universal gauge couplings gi = g in the kinetic terms ()5.1|) 
of the scalars one finds the mass terms for the gauge bosons 

N 

g 2 [v 2 A iti A? + u 2 (A t , - A {l+l) ,f] . (5.3) 
i=i 

After diagonalisation, the mass eigenvalues M n of the gauge bosons read 

M 2 n = g 2 v 2 + 2g 2 u 2 (l - cos 2vr^) , n = l,...,N, (5.4) 

which can be interpreted for n < Af as a linear KK spectrum ~ n/R = gu ■ n/N with a 
mass scale u ~ 10~ 2 eV. The scalar fields $j provide in addition a constant bulk (or kink) 
mass of the order v ~ 10 2 GeV. At this stage we can already observe the correspondence 
with the discretised boson spectrum from Eq. (|4.2|) . by identifying the VEV u with the 
inverse lattice spacing a -1 , and respectively the VEV v with the bulk mass M 8 . 

Furthermore, we find in £ mass [\l/j, Qi] from Eq. ()5.2|) terms of the type Qj^I/jL^fi+im + h.c. 
which yield fermion masses of order u when Qi has taken on its VEV. In the next section 
we will see that the ^i mass spectrum corresponds to the fermionic mass spectrum found 
in Eq. (|4.17|) . where the fermionic Casimir effect was calculated. Also in this case we will 
be able to identify u with the inverse lattice spacing a -1 in the ED. 

It is also possible to interpret the set of scalars $j as one 5D massive scalar on a trans- 
verse lattice, since in the potential V there exist a bulk mass term and terms of the 
type Y^i=i M 2 \&i — Qi&i + i/u\ 2 that can be interpreted as the discretised version of the 
continuum KK mass term j Q dy(d^) 2 . The resulting KK mass spectrum corresponds 
to the one for the gauge bosons in Eq. ([5.4)1 . but with an inverse lattice spacing of the 
order Mb and a constant kink mass Mq> ~ 10 2 GeV. 
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Let us close this section with a summary of the relevant scales that follow from this 
model. For the gauge fields A ifJi and the fermions we find the inverse lattice spacing to 
be of the order u ~ 10 -2 eV with an additional bulk gauge mass term of the order v ~ 
10 2 GeV. The fermions do not have a bulk mass term. Finally, the scalars with a bulk 
mass Mj, ~ 10 2 GeV yield a lattice spacing given by Mb ~ (10 4 . . . 10 5 ) eV. 



5.2 Vacuum Energy in Deconstruction 

In this section we show that by applying the correspondence between gauge theories in ge- 
ometric and in deconstructed higher dimensions, it is possible to transfer the methods for 
calculating finite Casimir energy densities in higher dimensions to the 4D deconstruction 
setup. One therefore obtains an unambiguous and well-defined prescription to determine 
finite vacuum energies of 4D quantum fields which have a higher-dimensional correspon- 
dence. We will demonstrate this procedure explicitly with the deconstruction model given 
in the previous section, which finally yields a 4D vacuum energy density that is comparable 
with the observed value p t, s ~ 10~ 47 GeV 4 . 

Let us start with the energy density p of the 4D quantum fields in deconstruction with N 
KK modes which is schematically given by 

N r i 

pcx^ / dV y/fP + M*, (5.5) 

n=l J 

where the masses M n depend on N, R, and the spin of the fields. Without any knowledge 
of the fifth dimension, it would not be clear how to put these UV divergent expressions 
into a sensible (finite) form. However, since we now interpret the KK tower in terms of an 
underlying higher- dimensional theory with certain boundary conditions, the 5D Casimir 
effect provides a well-known procedure to handle these UV divergences in four dimensions 
and yields a finite result. 

We can make the correspondence even more definite by comparing the zero-point modes in 
Eq. ()5. 5j) . which follow from the deconstruction mass spectra, with the energy momentum 
relations (j4.2j) and (j4.17|) appearing in the Casmir calculations. In the previous section we 
have already found this relation explicitly by confronting the gauge boson/scalar masses 
from Eq. (jOj) . 

M 2 =g 2 v 2 + 2g 2 u 2 (l - cos 2vr^) , 

with Eq. ijOl . 

m 2 = 2a- 2 (l - cos2vr— ) + M 2 

N 

Here, we directly observe that the inverse lattice spacing a -1 corresponds to gu, and 
respectively the bulk mass M s to gv. For simplicity, we treat all bosons like real scalars 
since they differ only by their number of degrees of freedom, and for our purposes even 
this small factor 2 or 3 can be neglected. In this case the overall prefactor in Eq. ()5.5|) 
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is exactly given by that of Eq. (|4.5J) after integrating out the ED. This shows that the 
formal correspondence is exact for the bosons. 

Now, we show this correspondence also for the fermions Let us first mention that in 
this case the prefactor in Eq. (|5.5|) also exactly matches that of Eq. (j4.18|) after integrating 
out the ED. To determine the mass spectrum M n for the Dirac spinor = (^f n h, ^ n R) T 
we start with Eq. (I5.2jl . where the link fields Qi have acquired their VEVs u, 



N 



C mass [%, Qi] -> - • ^nL [* (n+l)R ~ *(«-1)r] + h.C. 



n=l 



This sum contains the boundary terms — f ^il^or and f ^/vl^(at+i)r, and \1/or and ^(jv+rjR 
are defined by 

^+1 = 7^.*! and q = T m -q mN , m G Z, 

where we distinguish between untwisted (T = +1) and twisted (T = — 1) fermionic fields. 
Note that this is the discretised version of the continuum boundary condition 

V(y + mR) = T m ■ V(y), m G Z. 

Now, the Lagrangian in matrix form reads 

N 



C m3SB [^i] = ^nLM nk ^ kR + h.C. 



n,k=l 



where the mass matrix M nk and its square M 2 = MM^ = M^M are explicitly given by 



M= U - 
2 



+1 
-1 +1 



-T 



+T 



-1 
-1 



+1 




, M 2 = T 



2 


-1 



— T 




0-1 -TO 
2 0-1 — T 

2 -1 



— T 



■10 2 
-1 
-1 



-1 
2 
2 



The squared masses m 2 of the fermions are found to be the eigenvalues of M 2 . Thus, the 
mass spectrum for untwisted fields reads 



u sin 27r- 



n 

TV' 



n = 1, . . . , N, 



and for the twisted fields we obtain 



u 2 sin 2 27r 



n 



N 



n 



,N, 



(5.6) 



(5.7) 



which is consistent with the spectra found in Ref. [63]. Note that only for odd N, both 
spectra become identical, which has already been discussed in the context of the odd-even 
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artefact in Sec. 14.41 With these mass spectra the Casimir energy of the fermions is (—8) 
times the Casimir energy of a real scalar, whereas in a non-lattice calculation the energies 
differ only by a factor of (—4) (see Sec. 14. 4j) . This comes again from the well known 
phenomenon of fermion doubling in lattice theory. 

Note that the fermionic spectrum given in Eq. ()5.6|) (and respectively Eq. (|5.7J) for twisted 
fields) is exactly that of Eq. (|4.17|) . 

m 2 = oT 2 sin 2 2vr— + M f 2 , 

where the inverse lattice spacing a -1 is identified with the VEV u and vanishing bulk 
mass Mf. 

We therefore conclude that with the prescription motivated in this section, we are able 
to assign the finite Casimir energy density as calculated in Sees. 14.21 and 14.31 to the vac- 
uum energy densities of 4D quantum fields in deconstruction. Before we discuss the 
phenomeno logical consequences of this result, let us first treat the problem of fermion 
doubling. 

One way to handle the fermion doubling is to add a Wilson term [72] to C mass \^i,Qi\. 
This leads to a modified fermion mass Lagrangian, 



N r 



U 



£ 

n=l 



777 (<$ 

U 



^nL ( — #(n+l)R - * 



(n-l)L 



(5. 



which yields with (Q n ) = u the mass matrix M. The squared masses m 2 of the fermions 
are now given by the eigenvalues of M 2 , which is exactly the mass matrix for the bosons: 



M = u ■ 



-1 



-1 



+T 







M 



u 



2 

-1 



-T 



-1 

2 



-1 
-1 



2 

-1 



—T 



For the untwisted fields we get 



2u 



-2 



n 

COS27T-— 



n 



1 ... N, 



(5.9) 



and for the twisted ones n is replaced by n — |. These mass spectra are identical with 
the mass spectra of real scalars and yields the usual (continuum) factor (—4) in the 
vacuum energy density between the two field species. Looking at the above mass spectra, 
we remark that the spectrum (|5.9J) for scalars and Wilson-modified fermions does not 
contain a zero mode in the case of twisted fields. 

Let us now discuss the overall value of vacuum energy density p originating from the 
deconstruction model of Sec. 15. 11 According to the discussion in Sec. 14.41 the gauge bosons 
and scalars give a negative contribution to the CC. Without bulk masses its magnitude 
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field 



scalars $, 



gauge bosons Ai 



Table 5. 



fermions ^, 
The 4D fields 



inverse lattice spacing a 



-l 



M b ~ 10 4 - 10 & eV 



u 



1(T 2 eV 



it ~ 10" 2 eV 
lie deconsli 



bulk mass M 



m 



10 2 GeV 



10 2 GeV 







MR 



N ■ 10 b 



N ■ 10 



13 







in Llie UeconsLi ucLion scenanu ul 5ec. IO.ll wlncli can be hiLerpreted as 
KK modes of a 5D field on a transverse lattice. The bulk mass is denoted by M, 
and R = Na is the circumference of <SjL with N sites and a lattice spacing a. For 
large values of MR the Casimir effect will be highly suppressed according to Fig. 14.31 
and Table l4~21 Only the fermions contribute significantly to the CC. 



would of the order R~ A , where the size R of the ED is given by R = Na. For a small 
number iV of lattice sites this would induce a large negative contribution \p\ ^> p Q b s 
to the CC due to the small lattice spacing a -1 = Mb ~ 10 4 " 5 eV for the scalars $j. 
Fortunately, the scalars and also the gauge bosons are equipped with large bulk masses, 
which sufficiently suppress the Casimir energy thereby avoiding serious problems with 
observations. In contrast to this, the fermionic fields with KK masses of the order of 
the small VEV u ~ 10~ 2 eV induce a positive contribution {u/N) A ~ (10~ 3 eV) 4 to the 
CC which is of the observed order of magnitude already for a small number iV = 0(1) of 
lattice sites. In TableEHJwe give the relevant quantities for the vacuum energy suppression 
of the fields. Finally, it should be noted that we have determined only the vacuum 
energy contributions of quantum fields in the deconstruction setup which have a higher- 
dimensional correspondence. Other sources of vacuum energy could still lead to a CC 
that is in contrast with its observed value. Further discussions about dark energy in the 
context of discretised space-times and theory spaces can be in Refs. jT^l EH ES] ■ 



6 Dark Energy and Discretised Gravity 



Up to now, we have ignored gravity in the ED since in the usual approach of decon- 
struction gravity is assumed to have completely decoupled from the matter fields. In a 
gravitational ED, however, the gravity field has a discrete structure, too. By transferring 
the methods of deconstruction to this new setup one obtains an effective theory of a dis- 
cretised gravitational ED [HUES]- One has found that the gravitons that arise in such a 
framework exhibit a non-trivial feature in the form of a strong coupling behaviour that 
might occur significantly below the usual scale, where the effective theory breaks down. 
In addition, this new strong coupling scale depends on the size of the ED thereby leading 
to a so-called UV/IR- connect ion. In this chapter we will investigate the vacuum energy 
arising from the Casimir effect of massive quantum fields, where the strong coupling scale 
sets an upper limit for the bulk field masses. To avoid Casimir energy densities above the 
observed value of the CC large bulk field masses are needed for small EDs, but here they 
cannot be arbitrarily large due to the strong coupling. Because of this feature and the 
UV/IR-connection we are able to derive a lower bound on the size of the ED |76j . 

In the next section we will briefly introduce the effective theory of one discrete gravita- 
tional ED and the strong coupling scale arising in this model. Afterwards we will discuss 
the Casimir energy in this setup and derive the bounds on the fifth dimension. 



6.1 Gravitational Extra Dimensions 



Since we will perform the discretisation of a 6D gravitational model in great detail in 
Chap.[7|we keep the discussion of the 5D case relatively short. In the following we therefore 
present mainly results from the literature. For the 5D case of this chapter, an effective 
theory for a single discrete ED has been recently proposed |8T|l32] by implementing gravity 
in theory space similar to the concept of deconstruction. In this model gravity remains 
continuous in the bulk. 

Let us consider the theory space as given by the moose diagram shown in Fig. 16.11 Each of 





Figure 6.1: Part of the gravitational theory space for a discrete nfth dimension compactified 
on the circle S . Each site corresponds to one general coordinate invariance GCi 
(i = 1, 2, . . . , N), where two neighbouring sites i and i + 1 are connected by one link 
field Yi and we identify i + N = i. 
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the sites shown as circles corresponds to one general coordinate invariance (GC) symmetry 
GCi with % = 1 . . . N and can be identified with the 4D metric g % for this site. Every 
site i is connected to its neighbouring site i + 1 by the link field Y ifl , which transforms as a 
vector under the two neighbouring GCs. This is shown as an arrow connecting the GCs in 
the diagram. Moreover, the theory space and thus the ED is compactified on a circle by 
the identification i + N = i. Finally, we have on each site the usual 4D Einstein-Hilbert 
action as described by 




S 9 Bite = 2^ M 4 / d^vWW), (6-1) 

i=l 

where Rm{g l ) is the Ricci scalar on the site i, and M4 denotes the universal Planck scale on 
the sites, which is related to the continuous Planck scale Mpi = l/y/8irG by M| = MpjN 
with G as the 4D Newton's constant. From Eq. (j6.1|) we observe, that the action 5f ite 
is invariant under the product group Hf =1 GCi, which is explicitly broken to the diagonal 
GC by the gravitational interactions Sf ink between the sites. In our minimal discretisation 
with only nearest neighbour interactions, the action 5^ nk has the form a Fierz-Pauli 
mass term 1 

N 

SLk = E M " /d 4 *^™ 2 ^-^ 1 )^ ( 6 - 2 ) 
i=i ^ 

where the mass m corresponds to the inverse lattice spacing a -1 . The size of the fifth 
dimension is therefore given by R = N/m and the 5D Planck scale by M 5 = (M^/i?) 1 / 3 , 
which defines the usual UV cutoff of the 5D theory. Let us now expand in the weak 
field limit the metrics g* around flat space, g l = 77^ + /i' , where rj^ is the Minkowski 
metric. This leads to Fierz-Pauli graviton mass terms given by 

Sff = J d 4 x M 2 m 2 (25 M - 6 iJ+1 - ~ W)- ( 6 - 3 ) 

Note that this is the same mass matrix as for the gauge bosons and scalars in Sec. 15.11 
implying the graviton mass spectrum 

Tl TX7X 

m 2 n = 2m 2 (l - cos(2vr— )) = 4m 2 sin 2 (— ), n = 1, 2, . . . , N. (6.4) 

This spectrum describes one diagonal zero-mode graviton which corresponds to the unbro- 
ken GC and a phonon-like spectrum of massive gravitons that approximates for n <ti N 
a linear KK tower. At this level, the phenomenology of the model appears to be very 
similar to that of a deconstructed gauge theory. An important qualitative difference to 
deconstruct ion, however, reveals itself in the peculiar strong coupling effects of the theory. 

In Ref. jST] it was shown that the strong coupling behaviour is most conveniently dis- 
cussed by making use of the Callan-Coleman-Wess-Zumino formalism for effective field 

lr The Fierz-Pauli form for graviton mass terms ensures the absence of ghosts in the spectrum. For a 
recent discussion of ghosts in massive gravity, see Ref. [7S]. For other aspects in this context see, e.g., 

Refs. [23 HSU- 
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theories [EI]. Following this lead, the product symmetry group Tlf =1 GCi, which is broken 
by Sf ink , can be formally restored in Sf ink by including Goldstone bosons. Therefore, one 
expands each link field around x M as = x M + 7rf , where the Goldstone bosons trans- 
form non-linearly under GCi and GCi + \. The massive Goldstone vector bosons, which 
have three degrees of freedom, are eaten by the massless gravitons with two degrees of 
freedom, to generate the five polarisations of the massive gravitons with the spectrum 
given by Eq. ()6.4|) . Including the Goldstone bosons, the link action therefore acquires the 
form 



N 



(6.5) 



x (i-WWl^-W^i 1 ). (6.6) 



i=l 

u uJ-i y 7 5 ; \y a p u a^i ^pj-i y pa 

Let us further decompose the Goldstone bosons 7rf into its transverse and longitudi- 
nal (fii components as 7rf = Af + <9 M 0i. By using this form in Eq. ()6.5|) one finally obtains 
the interaction terms for the 0j that are responsible for the strong coupling, they have the 
form 

SL k = ...(dW0)(<9 2 0). (6.7) 

Actually, the interactions of the lowest lying scalar longitudinal component of the Gold- 
stone bosons lead to scattering amplitudes that quickly grow with the energy E of the 
scalars and lead to unitarity violation once the strong coupling scale is reached. For the 
present model it was found that the amplitude A(4>4> ~~ * 00) f° r — scattering is of the 
order A ~ E w /A\°, where 

M P1 N V 



is the strong coupling scale of the theory, that is set by the triple vertex of as described 
by Eq. ()6.7j) . From Eq. (|6.8J) . it is seen that the UV cutoff scale A 4 of the effective 
theory depends on the IR scale R of the compactified ED. This phenomenon has been 
called UV/IR connection [^21 since in a sensible effective theory for massive gravitons the 
lattice spacing a = m~ l must always be larger than the minimal lattice spacing defined 
by a m in ~ . This implies that the theory does not possess a naive continuum limit. 
In other words, for a given radius R, the effective theory is characterised by a highest 
possible number of lattice sites iV max = i?A 4 , which limits how fine grained the lattice can 
be made. 

In addition to the triple derivative coupling of 0, the Goldstone boson action contains other 
types of vertices, each of which can be associated with a characteristic strong coupling 
scale for that interaction jH2]. As two such typical examples, we will consider the scales 

Mpi\ 1/3 , /M P1 xl/5 



which we will later compare with A4. It is important to note that here the existence of the 
strong coupling is qualitatively different from the UV cutoff in deconstructed gauge the- 
ories. The strong coupling scale in deconstruction associated with the non-linear sigma 
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model approximation is of the order the inverse lattice spacing. Therefore deconstruc- 
tion may provide, unlike the effective theory of massive gravitons discussed here, an UV 
completion of higher- dimensional gauge theories. 



6.2 Bounds on the Size of the Extra Dimension 

Let us now investigate the Casimir energies of matter fields propagating in the discrete fifth 
dimension introduced in the previous section. For this purpose, we treat the gravitational 
theory space as a flat background for quantum fields propagating in the space-time with 
a discretised fifth dimension. Since these Casimir energy densities contribute to the CC, 
we require that they lie below the observed value p Q b s ~ 1CT 47 GeV 4 , associated with the 
accelerated expansion of the universe. In Sec. 14.21 we have found that for massless bulk 
fields the 4D Casimir energy density p scales with the size R of the ED as \p\ ~ R~ A , 
which would lead to a lower bound 2 R > (10 _3 eV) _1 ~ 0, 1 mm. A much smaller size R 
becomes possible, if the bulk fields have nonzero masses M, which implies according to 
Sec. 14.51 an exponentially suppression for M ^> R" 1 . In the discrete gravitational EDs, 
this suppression is only limited by the strong coupling scale A of the theory, since in a 
sensible effective field theory M should be smaller than the UV cutoff A. By virtue of the 
UV/IR connection, however, the cutoff A depends on R and can be much lower than the 
usual 4D Planck scale Mpi ~ 10 19 GeV. As a consequence, we expect from the Casimir 
effect a smallest possible value or lower limit on the size R, when M can at most be as 
large as the strong coupling scale A. 

In the ED, the boundary conditions for the quantum fields can be periodic or anti-periodic 
corresponding to untwisted and twisted fields, respectively. We have found in Sees. 14.21 
and 14.31 that the Casimir energy densities of these field configurations differ by a small 
factor and have opposite sign (ignoring the odd-even artefact). Following Eq. (|4.12|) . the 
4D Casimir energy density of a single untwisted real scalar field in the discretised fifth 
dimension can be written as 



1 4?r 

Puntwisted 7T7r> \q ' ~~rT 
1(M) 6 8 



N i 

m 4 . In m n — N ■ I ds ■ m\ In m s 

n=i Jo 



(6.10) 



where the ED has been integrated out. The mass spectrum m n is given by Eq. ()4.2j) 
and reads rn? n = 4m 2 sin 2 (im / N) + M 2 . Furthermore, the variable s is treated in the 
integral as a continuous parameter which replaces n/N in the sine function. As long as 
the number of lattice sites is N > O(10), the Casimir energy density on the transverse 
lattice in Eq. ()6.10|) differs less than < 1% from the value in the naive continuum limit 
N — > oo, which is demonstrated in Fig. 14.31 In the remainder of this section, we will 
therefore employ the expressions for the Casimir energy densities of quantum fields in 



2 A scenario for obtaining the observed CC from a 5D Casimir effect of massless bulk matter fields with a 
sub-mm extra dimension has been proposed, e.g., in Ref. [HSj- Current Cavendish-type experiments, 
however, put already very stringent upper bounds of the order R < 0, 1 mm on the possible size R of 
extra dimensions |84) . 
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the continuum theory as discussed in Sec. 13.21 In this approximation, the vacuum energy 
density of a real (un) twisted scalar field reads 

±1 (2vr) 5 r° (n 2 - x 2 ) 2 

P(un)twi S ted - g (27r)2 RA dn exp{27Tn)±V (P.ll) 

where the plus and minus signs belong to twisted and untwisted fields, respectively, 
and x = MR/(2n), in which M denotes the bulk mass of the scalar field. The integral 
in Eq. (|6.11|) can be performed exactly after neglecting the term ±1 in the denominator. 
Hence, both densities differ only in an overall sign: 

{MRf + 3MR + 3 MR 

P(un)twisted = ± (2tt) 2 R^ ' ^ ' 

When taking the sum of contributions for twisted and untwisted fields, the integrals must 
be added before carrying out the approximation, which gives 

4(MRf + 6MR + 3 2MR 
Psum " 16(2tt) 2 ^ e • ^ A6 > 

The corresponding energy densities of Dirac fermions are obtained by simply multiplying 
the scalar densities p( U n)twistod by —4, where we assume that the fermion doubling has been 
taken care of, e.g., in the way discussed in Sec. 15.21 Note that the applied approximation 
works fine even in the limit of vanishing bulk masses M — > 0. The basic feature expressed 
in Eqs. (I6.12J) and ()6.13|) is that for large bulk masses M > R- 1 , the energy density of 
massive matter fields becomes exponentially suppressed, which compensates for the large 
factor R~ A when R is comparatively small. 

Now, we are in a position to calculate the Casimir energy densities with the bulk masses M 
set equal to the strong coupling scales A 3 , A 4 and A 5 given in Eqs. ()6.8|) and (j6.9j) . The 
effective field theory description suggests that these are the largest possible values that M 
can take in the gravitational theory space. If the UV cutoff A is much larger than ~ -R -1 , 
the expressions in Eqs. (j6.12|) and (j6.13|) are dominated by the exponential damping 
factors, such that the Casimir energy densities are most strongly suppressed when M 
becomes of the order the strong coupling scale A, with A = A3,A4,As. Moreover, from 
Fig. 14.31 and Sec. 14.51 we know that this suppression is most effective, when the number of 
lattice sites N is maximised. Therefore we choose the inverse lattice spacing m = N/R 
to be also of the order A. 

Now the lower limit -R m i n on the size R of the ED emerges from requiring that the Casimir 
energy densities remain below the observed value p Q b s ~ 1CT 47 GeV 4 of the dark energy 
density. The results for an untwisted scalar field and the sum of twisted and untwisted 
fields are plotted in Fig. 16.21 Since the smallest value -R m i n that R can take is due to the 
UV/IR connection a function of A, we have considered .R mm (A) for all three scales A = 
A3,A4,A5. These values together with the corresponding maximum number of lattice 
sites = -R m i n ■ A(i? min ), where A(i? min ) is the strong coupling scale associated with Rmin, 
are summarised in Tab. 16.11 Note that we can apply here the relations from the continuum 
theory, since the number A^ of lattice sites is of the order 10 2 . Furthermore, the lattice 
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untwisted 




A(-R m in) 


N = -R m in • A(i? min ) 


A 3 


6, 1 • 1(T 12 GeV" 1 


3,6- 10 13 GeV 


219 


A 4 


9,0 ■ 10- 1U GeV" 1 


2,2- 10 11 GeV 


198 


A 5 


1,1 ■ lO^GeV" 1 


1,7- 10 9 GeV 


179 



sum 



A, 



A, 



A, 



R 



mm 



2 • 10" 13 GeV 



6,6- 10 



4,4-10 



11 GeV' 1 
" 9 GeV" 1 



1,4- 10 14 GeV 



1,6 ■ 10 12 GeV 



2, 1 • 10 10 GeV 



N = R„ 



A(R T 



112 



103 



95 



Table 6.1: The lower bounds R mm on the size R uf Llie ED fur an untwisted real scalar field 
and the sum of a twisted and an untwisted scalar. Additionally, the values of the 
strong coupling scale A and the number of lattice sites N are given when R is equal 
to -Rmm- For the scale A, we considered each of the three choices A = A3, A4, A5 from 
Eqs. Ijfi.HJI and The lower bound i? m i n emerges from the requirement that the 

absolute Casimir energy density lies below the observed value p bs of the DE density, 
when the bulk field mass M takes the largest possible value M ~ A. 



calculation leads to energy densities (drawn as circles in Fig. I6.2|) . that agree very well 
with the values in the continuum theory. For -R m i n the values of the continuum and lattice 
formulas differ by about 15%, which is negligible, since the strong coupling scales A 3) 4 i5 
from Eqs. (|6.8|) and ()6.9|) are order of magnitude estimates. For instance, the lattice 
calculation for an untwisted scalar field and A = A 3 gives Rram — 6.8 • 10 _12 GeV _1 , 
whereas the continuum approximation yields R m i n = 6.1 • 10~ 12 GeV -1 . 

Let us now discuss the results. For the sum of a twisted and an untwisted field, we observe 
that the Casimir energy density of massive bulk fields exhibits a stronger suppression due 
to the different signs of both components. From Fig. 16. 2\ we read off that the minimal 
radius -R m i n of the discrete gravitational extra dimension lies in the range 

(KPGeV)- 1 < R min < (10 7 GeV)-\ (6.14) 

where we typically find A(i? min ) ~ 10 2 x R m ] n . For a radius R which is much smaller 
than the range given in Eq. ()6.14|) . the Casimir energy densities of the bulk matter fields 
would significantly exceed p Q b s and thus run into conflict with observation. Of course, 
there may be other possible sources of dark energy which might be responsible for the 
accelerated expansion of the universe, but it seems unlikely that they could exactly cancel 
the potentially large contributions from the Casimir effect in EDs. 
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R/GeV" 




R/GeV" 




Figure 6.2: For the three 6 choices A = A3, A4, A5 of the strong coupling scafe A from Eqs. (|6.8|) 



and (|6.9j) . we have plotted the values of A, the Casimir energy densities p, and the 
corresponding number N = RA of lattice sites as functions of the size R of the 
fifth dimension. The energy densities p are given for the untwisted scalar field, see 
Eq. (|6.12|) . and for the sum of one untwisted and one twisted scalar field as given 
by Eq. (|6.13|) . Note, that p is negative in both cases, and the bulk masses of the 
fields have their maximal values, given by A, according to Sec. 16.21 In the plots of p, 
the horizontal dashed line marks the observed value p \, s ~ 10 -47 GeV 4 of the DE 
density and the circles represent exact lattice values from Eq. (|6,1U|) , 



7 Discretised Curved Disk 



This chapter is devoted to a special 6D model, where as before the 4D subspace is a flat 
Minkowski space-time. Both higher dimensions, however, form the discrete version of a 
disk with constant curvature. In the case of positive curvature, for example, the disk would 
be part of a 2-sphere. Moreover, we apply a very special discretisation, where on the disk 
boundary we place N lattice sites and only one site in the centre of the disk (Fig. 17. 1)1 . 
For gauge theories and gravitational EDs similar setups have been investigated before in 
Refs. j66j|85], but here we include the possibility that the disk has a constant curvature, 
too [HE]- This more general case allows a flexible hierarchy in the mass structures of 
the effective 4D fields, which can be controlled entirely by the geometry parameters of 
the disk. Another motivation for the setup of this chapter was found in the context of 
multi-throat geometries, where the possibility of hiding large EDs was discussed j87j . 

In the following we will first discuss the 6D model in the continuum and show that it is 
a solution of Einstein's equations. Afterwards we explicitly derive the discretisation of 
the EDs and determine the mass spectrum for the effective 4D gravitons and the value 
for the effective 4D Planck mass. All these scales will be found to depend directly on the 
parameters of the curved disk. The explicit calculations also help to understand the 5D 
case of Sec. 16.1) especially the origin of the Fierz-Pauli graviton mass terms in Eq. ()6.3j) 
and the link action ()6.2j) . Finally, we show how to implement fermions into this scenario, 
where we will find that the corresponding mass spectrum for the 4D fermions is directly 
related to the graviton spectrum. With these results we finally show a simple application 
to generate small fermion masses. 

7.1 Curved Disk Geometry 

Let us consider an effective field theory for massive gravitons, which arise after com- 
pactification of two gravitational EDs that have been discretised. We start out with a 
continuum theory, where 6D general relativity is compactified to four dimensions on a 
two-dimensional disk of constant curvature. In this space, the 6D coordinates are denoted 
by x with capital Latin letters M = 0, 1, 2, 3, 5, 6, while Greek letters appear in the usual 
4D coordinates x^ as /x = 0, 1, 2, 3. The position of a point on the disk is described by 
the radial coordinate r := x 5 G [0, L\ and the polar coordinate (p := x 6 G [0, 2tt]. The 6D 
Minkowski metric is given by r] MN = diag(l, — 1, . . . , — 1). In the following, we consider a 
metric qmn that is defined by the line element 

ds 2 = g„J x^dx^dx" ^dr 2 - r 2 d<^ 2 , (7.1) 

1 — er 2 
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j=N=12 



j=io 




\ j=4 



Figure 7.1: The curved disk model of Chap. with 1 a special discretisation, where N lattice sites 
are located on the boundary and only one site (j = 0) sits in the centre. In this 
example we chose N = 12. In the effective theory of massive 4D gravitons the solid 
lines in radial direction correspond to the mass scale m* from Eq. (j7.13[) . whereas m 
characterises the interactions in angular direction, which are shown as dashed lines. 



where 1/ y|e| is the curvature radius of the disk and g fMU (x M ) the metric of the 4D sub- 
space. Note that the parameter e controls the curvature of the disk, for e > the disk 
is spherically curved, and e < leads to a hyperbolic disk 1 , respectively, e = corre- 
sponds to a flat disk. From Eq. (|7.1|) . we read off = — (1 — er 2 ) -1 and = — r 2 . 
We denote partial and covariant derivatives by commas and semicolons, respectively. Us- 
ing (755^ = — 2er/(l — er 2 ) 2 and (766,5 — ~ 2r we find for the non-zero Christoffel symbols 

4D 1 1 1 
-pa _ -pa pM _ _ n W n _ _ n W n p5 55 

1 iiu - ~2~9 fiV,6, J- 55 - 2~9 #55,5 - ^ _ 2 , 

^ 6 = -^W = -r(l-er 2 ), r« 6 = ^%,5 = p (7.2) 

where T^:= \g rjp {g i i P , v + — 9yv,p)- With our conventions the Riemann (curvature) 
tensor R A MND , the Ricci tensor Rmn and the Ricci scalar R are defined by 

pA _ -pA -pA i pA -pB -pB -pA 

n MND ■— 1 MD,N 1 MN,D 1 BN MD 1 MN BDi 

p pA -pA _ -pA , pA p_B _ p_B pA 

" MN •— '> MN A — 1 MA,N 1 MN,A "r 1 BJV AfA 1 MN BAi 

R '■= Rmn9 MN ■ 



1 Compact hyperbolic extra dimensions have been discussed, e.g., in Refs. |88l I89j . 
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The 6D Einstein-Hilbert action 2 including a CC A and some matter action Setter reads 



S = M A & J d 6 x y/\g[(R - 2A) + 5 matter , 



(7.3) 



where g = det Qmn an d M 6 denotes the 6D Planck scale. Note that in this chapter we 
denote the CC by the quantity A, which has mass dimension two and is related to the 
vacuum energy density A via the Planck mass. The requirement that the variation SS 
with respect to the metric vanishes, leads to Einstein's equations 



Rain — ^9mnR + ^9mn 



2M 6 4 



T 



MN, 



and the energy- momentum tensor Tmn follows from 

(5 Scatter 



MN 



Sg 



MN 



Now we are interested in the Einstein-Hilbert action corresponding to the specific 6D met- 
ric given by Eq. (|7.1)l . Under the 6D integral we therefore calculate \J\g~\R = \/\g\(g^ P R fJ ,u+ 
g 55 Rn5 + g m R§§) in the following. The first term reads 



nln^R - * An\n liU IT 1 " — V a 4- V a T@ — T a 
y\y lh nv — V \y\y l fJ.a,u 1 tiu,a ' 1 /3i/ L fia 1 fiv 1 f3a 

p5 I p5 p/3 , -pa p5 p5 pA 

1 fj,v,5 "r" ^5 i" 1 5V 1 fia 1 /xi/ 1 5A 

p6 I p6 p/3 , pa p6 p6 pA ] 

1 1 /g^ 1 ^6 ' 1 6V 1 /ita 1 fj,u L 6Aj • 

The derivative of the Christoffel symbol with respect to r (and similarly with respect to tp) 
can be written as 



J ,5 



where on the right-hand side we have used (yig|),A = vM^ba i n the second term 
and g IJ,u 5 = —g m g v ^g a ^^ = —2g^ a T u ab in the third term. Thus we find 



\g\g" u R 



fJbV 



\g\R. 



ID 



g\g^Tl 



\g\!TT% 



2 In n space-time dimensions the Einstein-Hilbert action with Planck scale M n reads 

S = J d n x^\Adi n - 2 \R - 2A) + S matter , 

yielding Einstein's equations in the form 

Rmn — ~9mnR + ^9mn = 7z 9 n Tmn ■ 

2 2MA ' 
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with i?4D denoting the usual 4D curvature scalar. Similarly we find for the other terms 



r>-> 



W b [r 



A 

5A,5 



\g\srr 



A5 



J ,5 



^A , -pA -pB 
55,A "r 1 B5 1 5A 



r B r A 1 

1 55 1 BA\ 



+ 



[21*1^ - rf 5 rl + r| A r 

op5 p5 I pA pB -pB -pA ] 

Zi 55 1 55 B5 1 A5 1 55 1 BA\ • 



55-pA 
A5 



g\g 00 r 



\g\9 55 r£ 5 



.A 



+ 



and 



I I 55r-pa p/3 -pa p/3 _ npa p6 l 



\ n \n m \r A -r A +r A r B -r B r A 1 

|y|y L 1 6A,6 L 66,A^ B6 L 6A 1 66 1 BA\ 



9\9 1 A6 



|y|y 1 66 



.A 



+ 



i pA -pB -pB -pA 1 

"i" 1 B6 1 A6 1 66 1 BAj • 



P_B -pA , pi? pA 
J- Rfi-L 4fi ~r 1 p a J- g 



BA- 1 - 66 



\g\g m r 



66-pA 
A6 



\g\g 



66 -pA 
66 



+ 



I I 66rpa p/3 _ pa p/3 1 

|y|y l 1 /36 1 «6 ^e 1 ^.! 



Within these terms we collect all total derivatives that read altogether 



\g\(-g^r% + g 6& T 



a6) 



+ 



g\( g'-vL ■ rr;,:, -- <nt, ■ g^l.) 



where we now write the last two terms as 



55"po p6 
a5 L 65 



Here, we have used the relations g m T^ 



-2^\g\g 55 rl 5 rt 5 . 

-g 55 Tl 5 and I* 



(7.4) 

66 - -y -65 — -65,5 - -( r 6 5 ) 2 following from 
Eqs. (17. 2j) . Therefore, the first term in Eq. ()7.4j) cancels one term in g 55 R^ given above, 
and the last terms yields the cosmological term y/\g\(2e). 

Finally, we are able to divide the 6D action ()7.3j) without matter into three parts, 



S = S, 



4D 



s. 



surface 



+ s n 



(7.5) 



where the (modified) Einstein-Hilbert action S^n for the 4D subspace, the remaining 
surface terms S^face an d some terms S ma , ss , that will later become graviton mass terms, 
are respectively given by 



Sad 

^surface 



Ml 



Mi 



Mi 



+ Mi 



dV 
d 6 x 
d 6 x- 
dV 



\g\{R 4B + 2e), 



I 55 pa 

y|y 1 a5 



+ 2 



glg m T * 



06 



I I rt 55rpa p/3 _ pa p 
\g\g l L /35 1 a5 1 05 1 



J 1 
ah 1 - /35J 



I ^Tpa p/3 _ po p 
y|y l L /36 1 a6 1 06 1 



/36J 



(7.6) 
(7.7) 

(7.8) 
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Let us first consider S , sur face, which reads in our setup 



surface 



2 77 



M 6 4 / d 4 x / d<p dr 



<?I<7 5 W,5 



+ 



J ,6 



The (^-integral over the second term vanishes if we apply periodic boundary conditions, 
which is plausible for a disk, and the remaining term 



^surface = M$ / d*X / d(f 



2tt 



r=0 



(7.9) 



can be also be removed by choosing appropriate boundary conditions at r = and r = L. 
In addition, the lower limit would vanishes because of \/\g\g 55 oc r as long as does 
not diverge at r = 0. 



Furthermore, we write the mass term action in a simplified form given by 



<7 - M 4 

'-'mass — 1V1 Q 



d 6 x^/\i\ f - ~ A 9 cc g,uM u 9 aP - g m g v0 )g 



c=5,6 



a/3,c 



(7.10) 



Note that the naively discretised version of this equation would become in the 5D setup 
of Sec. 16.11 similar to Eq. ()6.2|) . which will become more obvious in Sec. 17.51 



7.2 Solving Einstein's Equations 

In this section we will find a simple solution of Einstein's equation for the 6D metric given 
by Eq. (|7.1jl . For this purpose let us write the gravitational part of the 6D action ()7.3|) in 
the form 

S = Mt J d 6 x y/\g~\ [R - 2A • (g AB n AB )} , 

where we have introduced the tensor 

n AB := diag(0, 0, 0, 0, -g 55 , -g m ), 

given in the coordinates of the previous section. The term involving A is a modified 
cosmological term in the following sense. As long as we are working only with the action S 
as a number we can evaluate the expression g AB riAB = = 1 and thus obtain standard 
6D gravity with a 6D cosmological constant A. However, if we want to derive Einstein's 
equations from the action functional S^ab] we have to leave g AB riAB unevaluated. The 
variational principle SS/Sg AB = then gives 6D Einstein's equations 

Gain + 57*^ — Rmn ~ -^9mnR + g^u^ = 
with a CC A, that is located only within the whole 4D subspace described by the metric 



9fu, ■= 9mn ~ 2^a/jv = diag(sf o, 911,922, #33, 0, 0). 
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To calculate this we have used the relations 



Sty\g\R) 

Sg MN 



9\Gmn and respectively 



6(-2X 



g\g AB n A n B ) 



1 



MN 



-2X(-^y/\g\g M Ng AB n A B + V\d\ S M S N n AB) 



A zB„ 



— X(gMN — 2hmn) — svA. 

We now apply the simple ansatz g^ v = rj^ = diag(l, —1,-1,-1) leading to 

(-e + Xjrjftv 



Gmn + sy^A 











which becomes a solution of Einstein's equations GuN + g^v^ = by setting A = e. In this 
case A explicitly cancels the cosmological term in Eq. (|7.6)1 . In the next section we will 
investigate tensor fluctuations (gravitons) around the 4D Minkowski metric rj^, which we 
have found to be a solution of Einstein's equations. 



7.3 Massive 4D Gravitons 



In Sec. 17.11 we have found the terms ()7.10j) that will become, as we show now, the mass 
terms for 4D gravitons. Let us first write oVass in the form 



On 



1 



M 6 4 / d 4 x / d^ dr [ + -rVT^d r g^(g^g a P - g^g^)d r g aP 



+ 



a/3 



where we have used \J \g\ = \J \g±\ ^/g^g^ in the second line. Then we introduce graviton 
fields h^ u on a flat Minkowski metric r]^ u for the 4D subspace by the replacement 

g^v * hfiv- 

Note that in our approach we ignore graviphoton and radion excitations, which could result 
from the g 5 M and g S M components of the metric. Since r]^ u is constant we have g^A —> 
and by expanding S mass in second order in h^ u we find 



Ml J d 4 x / dipdr y + -rVl - er 2 d r h^rTi* ~ V^V^)d r h a p 



+ 



4 ryl — er^ 



To obtain the discrete disk as described at the beginning of this chapter we perform the 
discretisation of the disk in the following way (Fig. 17. 1)1 . We put iV points on the boundary 
and one point in centre of the disk so that only two points are lying in radial direction. The 
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coordinate distance between the centre point and every point on the boundary is given 
by the coordinate radius L of the disk, which in general differs from its proper radius. 
On the boundary we denote the graviton fields by h l with i = 1 . . . N, and the position 
is given by tp 1 — i ■ Atp, where Atp = 27i/N is the angular lattice spacing. The graviton 
field in the centre carries the index 0, and the lattice spacing in radial direction is 
given by Ar = L. Explicitly, we discretise the disk by applying the following replacements 

tO *SVna.ss* 



d r h(l{?) 

dr f{r) 



I 



(h* - h°) 

Ar" ' 

(h i+l - ti) 

£>r-/(L) = Ar-/(L), 

r=L 
N 



(7.H) 



i=i 



Note that the integral J dr is replaced by just one summation region of length L, where 
the summand is evaluated at the position r = L to avoid problems with the derivative d v 
at r = 0. Thus we obtain for the mass terms on the discretised disk 



5„ 



x 



Ml 



N 

d*x ^2 ^(pAr 
i=i 



+ -Wl -eL 2 ■ V V ( ^ a0 _ rfrfP) 
4 Ar i i j 



h 



a/3 



1 



1 



lA+1 



h 



[IV 



4 LVl - eL 2 



A<p 



Ar 

h l +} - h 



l a/3 



a/3 



Ay? 



which become more clear in the form 

N 



x 



M, 



J i=l 



(7.12) 



From the last line we observe immediately the relation to Eq. ()6.2j) that corresponds 
to the 5D case. Note that the actual graviton mass scale from the radial links, m*, and 
respectively from the angular links, m, depend on the 4D Planck mass M4 of the observer's 
site (brane): 



Mi 1 2tt 



M 2 4 



iV 



• Vl - eL 2 , 



m 



Ml 1 



N 



M\ 4 2vr VI -eL 2 ' 



(7.13) 
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However, the ratio of masses is independent of the Planck scales, 

m 2 (27r) 2 0n 

—2 = -7^- X - eL ' 7 - 14 

which also shows that arbitrarily large hierarchies between and m are possible by 
choosing eL 2 and N appropriately. 

Therefore the 6D action of the curved disk with A = e, 

S m = Ml J d e x VW\[R- 2K{g AB n AB )] , 
becomes with the discretisation procedure from above and the expansion — > r]^ + 

S 6D -> Mq J d 6 XA/|^|-R 4 D + Smass, 

where .Smass is now given by Eq. (|7.12J) . 

7.4 4D Planck Scale on the Sites 

Let us now turn to the determination of the 4D Planck mass M4 on the sites. For this 
purpose we need to know the proper area of the curved disk that is given by 



[ Vl-eL 2 ), e>0 

r 2 * r r L r e 

A := / d<p drV|c/5 5 ^ 66 | = 2vr / dr = <^ vrL 2 , e = 

J ° J ° J ° [fi(y/T+W?-l), e<0. 

(7-15) 

For \e\L 2 C 1 we find that the area scales like A = ttL 2 + 0(eL A ) and for e < the 
area grows approximately linearly with L for large values of eL 2 : A m 2ttL/ y/\e\. In 
the case e > the coordinate range r 6 [0,1/ y/e\ covers one half of a 2-sphere with 
radius l/\/e, thereby r = 1/y/e corresponds to the equator. The other half of the sphere 
can be described formally by the same metric and the coordinate range r e [—1/^,0], 
where r — > with r < means approaching the antipode. Thus the whole sphere has the 
usual area 4ir/e. Moreover, on the first half of the sphere the relation between the proper 
length s and the radial coordinate r is given by 



dr'V\9^\= / dr' 



1 e>0 



^parcsinfi/er) 



—= arcsinfi/er) 



'0 Jo Vl - er 2 

which is equal to the length of the arc on the right-hand side of Fig. 17.21 

For completeness we remark that in the hyperbolic disk case (e < 0) the proper distance s 
from the centre in radial direction is given by 



dr' — 

VI — er 



e<0 



2 



arsinh(y / [e[r) 



r 



arsinh( y^jefr) 
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Figure 7.2: Illustration of the relation between the coordinate r and the corresponding proper 
length s for the spherically curved disk with e > 0. The first half of the sphere is 
described by the coordinate range r E [0, l/y/e\, which covers the lower part of the 
sphere. The length s of the arc is given by a/ 'yfe with a = arcsin(y / er). The upper 
part is covered by the range r G [—1/y/e, 0], where a = arcsin(y / er) +7r/2 and r < 0. 



whereas for the flat case, e = 0, we simply find s 



r. 



We now proceed to discretise the term Mg J d^-y/jgji^D from the modified 4D ac- 
tion (|7.6jl . Since the extra-dimensional disk has a constant curvature it seems well mo- 
tivated that also the 4D Planck scales M4 should be constant and universal on all sites. 
Thus the terms of all N + 1 sites have to take on the form 



site i 



MlY^ I d 4 x v^#4 
Since we cannot allow R$e> to depend on r and ip, we thus have 



(7.16) 



(74I-R4D 



site i 



94 1 R 



"ID 



for all sites i. 



Now putting these parts together one finds that Mg J d 6 x^/\g\R4v has to be discretised 
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as 



N 

Ml [ d^drvWeel / d 4 x ^\R iD -> M\ ^ ! d 4 x ^\R^ (7.17) 



proper area A 



N 



! d 4 xv^i? 4D (7.18) 
Mj{N + I) J d 4 x y/\gT\R*D- (7.19) 



By comparing the left-hand side with the right-hand side the universal 4D site Planck 
scale is fixed, 

Ml = (7.20) 

where the proper disk area A as given in Eq. ()7.15|) is equally divided upon all sites. How- 
ever, it is important to note that M4 is not the (reduced) Planck scale Mpi = 1/ (8irG) ~ 
10 18 GeV that couples gravity to 4D matter. But Mpi is determined by integrating out 
the EDs in the continuum, which means 



M, 



e / d x V\9\R 



4D 



M, 



pi 



d x y\g4\R. 



l4D- 



From this we find as in the discretised 5D case of Sec. 16.11 the relation 

M|j = Ml A = (N + 1)M|. 



(7.21) 



In this section we have desisted from naively discretising the area integral J dtpdr a/ | #55 #66 1 
in S4D since ^/g^gm ^ r vanishes on the centre site and consequently there would not be 
a gravitational action. On the other hand, this strategy could be useful when considering 
a model, where the centre site just represents an auxiliary construction without physical 
relevance. However, we will not pursue this idea here any further. 

Instead, we mention another possibility, where one could obtain a hierarchy between the 
centre and boundary sites. Let us consider for this purpose the spherically (e > 0) curved 
disk with L equal to its maximal value on both hemispheres of the sphere. In this case, 
we have a closed 2-sphere, where the centre sits on one point and all the boundary points 
sit on the antipode. This means effectively that there are only two points left. It might 
therefore appear "natural" to assign both points (centre and boundary) one half of the 
total area 2A, which leads to 




We thus see that the squared 4D Planck mass on one boundary site would just be one N' 
of the squared 4D Planck mass in the centre. 
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7.5 Graviton Mass Spectrum 

Up to now we have determined the mass term action for 4D gravitons that follows from the 
discrete 6D model under consideration. Let us now determine the actual mass spectrum. 
First we apply the graviton expansion = r/^ + h^ v of the 4D metric also in the 
4D curvature scalar R^, that occurs in the gravitational terms Mj Yl!i=o I d^x V^T-^d 
for all sites. This leads to the kinetic terms for each graviton and was found to have the 
form [13 ED] 



9\R4D = V\g^\\{& i h VK d IM h VK - Vhdph - 2h»h, + 2h^h) + 0(h%), (7.22) 



where the following relations have been used: 



h",,, h. 



A'' 



hT := r\ m rf^K^ 



1 fiu 



-AW 



3 ) 

fiuJ 



Up to second order in we find that the whole graviton action, following from Eqs. ()7.22j) 
and 1)7.12)) . is given by 



S, 



graviton 



N 

a — n J 



2h tfl hl + 2h tfi d ll h t 



i=0 

A' 



x 



»=i 
-ml 

2 



where the mass scales and m have been determined in Eq. (|7.13j) . 



(7.23) 



Next we have to diagonalise the graviton mass terms by a unitary transformation. As 
mentioned in Sec. 17.41 we choose the site Planck scales M4 to be universal and fixed 
as in Eq. 1)7.20}) . Then the squared (N + l)-dimensional mass matrix M 2 in the ba- 
sis (h®, h*, . . . , hp V ) has the following form 



M 



(N 
-1 

-1 



1 



-1 



-1 
1 



-1\ 



+ m 



J 










2 



1 







-1 



-1 



-1 

2/ 



(7.24) 



where the first matrix proportional to originates from the interactions in radial di- 
rection. The second matrix comes from the angular interactions and is very similar to 
the gauge boson mass matrix described by Eq. (|5.3|) . Fortunately, both matrices can be 
diagonalised simultaneously by a unitary transformation. When we denote the graviton 
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mass eigenstates by corresponding to the mass M n , we find the following relations 
for the eigenvectors: 



fJbV 



I - 

± V y 



H 



N 



^^[sin(27ri#)+coB(27ri^)] -h 

N 



1=1 

1 



y/N{N + 1) 



i=l 



where p = 1 . . . N — 1. The eigenvalues M n are respectively given by 



Ml 
Ml 



M 2 N 







tt\ + 4m sm — , 
(JV + l)mJ 



(7.25) 
(7.26) 
(7.27) 



(7.28) 
(7.29) 
(7.30) 



From these results we observe that the zero-mode H® v has a flat profile and is equally 
located on all sites, whereas the mode H^ u with squared mass (N + l)m1 
on the centre site with equal support on the boundary sites. The N — 1 



is peaked 



modes if^, 



with p — 1 . . . iV — 1 are located only on the boundary with a typical discrete KK mass 
spectrum like Eq. ()4.2|) that has been shifted by ml. In the limit m the masses of 

the states from Eq. (|7.26j) become degenerate, and for iV 3> 1 the mode H^ v becomes 
very heavy. 



Finally, we remark that a scenario related to ours has been discussed recently in the 
context of multi-throat geometries jHZI- It was shown that large EDs can be hidden in 
the sense that the occurrence of massive KK modes is shifted to energies much higher 
than the compactification scale of the ED, which helps evading limits on KK particles 
and other bounds [25 ^ 126 } l27j. This behaviour can be observed here for the modes if 
in the limit m+ ^> m, too. 



n>0 



7.6 Fermions on the Disk 



Let us now investigate the incorporation of Dirac fermions into the discretised disk model 
of Sec. 17.11 As in the graviton case we start with a 6D Dirac fermion \l/ in the continuum. 
Using the vielbein formalism j^I], the corresponding action S on the curved disk reads IT(" 



S 



d*x J\g 



(7.31) 



where we denote 6D Lorentz indices by A, B , . . . and general coordinate indices by M, N, . . . 
respectively. Moreover, we call G A the 6D Dirac matrices, and for the barred spinor ^ we 
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use the abbreviation \1/ = \1/^G . The vielbein components V A M (x N ) follow from the re- 
lation g MN = V^V^t]ab, which connects the Lorentz coordinate system with the general 
coordinate system. For the diagonal metric given by Eq. (|7.1jl . 



ds 2 = g,Jx M )dx' 1 dx' y ^dr 2 - r'dip 1 

1 — er z 



(7.32) 



we find Vjm = $m w hh the exceptions V^ = l = 1/^55] and V^ = g = a/|^ 66 |. To avoid 
confusion with indices we denote the inverse vielbein components by 



V, := V, 



M=5 
A=5 



55 I 



and Vfi := V, 



M=6 

6 •- v A=& 



9" 



respectively. On a curved space-time the covariant derivative V m = du + Tm for spinors 
contains in addition to the usual partial derivative du also the spin connection 



M 



^[G A , G ]VaV B n 



(7.33) 



To determine the form of the 6D 7- matrices let us first look at the 4D case, where 
the 7-matrices are given by 



7 



' 






a k ' 







, l k = 


-a k 



«7o7i7273 



-1 
+1 



where k — 1 ... 3 and o k denote the Pauli matrices. This set of matrices has the properties 

(7°) f =7°, (7 fe ) f = -7 fe , (7 5 ) f = 7 5 , {7 5 ,7n = 

with the anti-commutator {□,□}. In five dimensions the number of spinor components 
is still four and the corresponding 7-matrices are simply given by r° = 7 , T k = 7 fc 

and r 5 = ir> = -(r 5 )t. 

In six dimensions, however, the Dirac algebra is 8- dimensional. Here, we use the following 
set of 7-matrices 



G° 
G n 
G 6 



1 4 
1 4 

r°r n 
-r°r n 

r° 
-r° 



(G°)t, 

1 =-(G") t , n = l,2,3,5, 
(G 6 )t, 



which fulfil the Cifford algebra {G A ,G B } = 2r] AB ■ 1 8 . Furthermore, in 6D one can also 
define a chirality matrix G c by 



r — Lx t_x Lx l_r l_r l_r 



-1 4 
1 4 



Gl 



which satisfies G\ — \ and {G c , G A } = 0. 
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From the form of the vielbeins and the 7- matrices it follows that the terms Vbn-m that 
occur in the spin connection 3 Y M vanish except for V(b=6,n=5)-,6 = 1 and Vm=5,N=6);6 — 
— ryl — er 2 , which leads to 



-[G 5 , G 6 ]^! - er 2 = -iV 5 



7 5 
7 5 



1 6' 



while all other components T M ^ 6 vanish. Finally, we observe that in the action (|7.31|) the 
term involving r 6 , 



vanishes because of r 6 ^G 6 = ^TlG G 6 = ^G 6 T 6 . 

Let us now diagonalise the action by the substitution \I/ := G 6 $, which yields \I/ 
—&G 6 G° and respectively 



-r° 
r° 



[G°d + G k d k + V 5 G% + V 6 G 6 d 6 ] G 6 $ 



in the action ([7.31)1 . It turns out that all the products G A G 6 are block diagonal and if we 
decompose the eight-component spinor $ = ($ a , $&) T into two four-component spinors $ a , 
$6 the last line reads 



i X 
+ 





° ^ 


(o° 


1° J 


—27 s 











d + 



-r 








<9a 



V 5 d 5 + 



7 
1 
-1 



x 



with = $^ b 7°. From this one can read off that $ a corresponds to but with negative 
energy, therefore we will work only with in the following. If we now denote the left- 
and right-handed components of $6 by $l,r : = |(1 = F7 5 ) < £b, then the full action for $ b 
can be written in the form 



S 



d b x^\g 



- V5- ($L<9 5 $R - $R<9 5 $ L - <9 5 $L$R + <9 5 $R$ L ) 

- iVg - ($l9b$r + 1r~<9 6 $ l - <9 6 $ L $ R - <9 6 $ R $ L ) 

In the last line let us further transform both middle terms by an integration by parts and 
apply periodicity conditions with respect to the (^-coordinate: 

/*2"7r /*27r 

/ dtp ($^9 6 $ L - d 6 $ L $ R ) = [$r"$ l - $Z$ R ] ^ - / dtp (Wl - $l<9 6 $ r ) . 
Jo Jo 



3 Note that we use subscript indices for the spin connection Tm and superscript indices for the 5D 
7-matrices T A . 
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And for both middle terms in the second line we respectively obtain 

/ dry 10550661^5 (-$r# 5 $l - d 5 $ L <l> R ) = [-r$^$ L -r$7> R ] 
Jo 

L 



[ dr (-$ R $ L - $ L $ R ) 
Jo 

/ dr • y^\^\ (-<9 5 $ R $ L - $7> 5 $r) 

</ 



By using the same discretisation prescription 

- [ dr (-$^$ L - $L$r) -> I [$R$L + $L$R] r=L 
./0 

as given in Eq. (|7.11jl for the gravitons in Sec. 17. 3| we find that only the terms in the last 
line will survive in the end. 

Since 6D spinors have mass dimension | we have to rescale them in order to obtain usual 
4D spinors. As in the graviton case we will integrate the kinetic terms over the EDs and 
also apply a similar discretisation procedure as in Sec. 17.41 which means 

3 

b- 



J - j=0 + J 



where A is the proper area of the EDs as given in Eq. (j7.15|) . Finally we absorb the 
factor A/(N + 1) into the fermion fields \ := ^y/Af (N + 1) and subsequently apply the 
discretisations 

d 5 X - (x J -X°)/Ar, 

d 6X -> (X j+1 -X j )/&P, (7-34) 

N 



! drdip-f - ^ArA^-/ 

J .7=1 



where we also use Ar = L and Atp = 2ir/N. Note that in radial direction we cannot use 
a symmetric derivative as in Sec. 14.31 since there are only two lattice points available. For 
the derivative in radial direction we will later on consider a more symmetric form. As a 
result we obtain the action for N + 1 4D fermions, 



s = J2 [ d4 4* (x j Yd„x j - d»x j Yx j ) 

j=o J 

N 

- / d * x ■ m * (%l(Xr - Xr) + (xi ~ Xk)xI / 
i=i J 

N » 

- ^2jd 4 x-i-m (xl{x£ l ~ xi) ~ Ur 1 - Xr)X 3 l) 



86 



Chapter 7. Discretised Curved Disk 



where the mass scales read 



m 



ArAipy/\g 55 g m \V* 



ArA(py/\g 55 g m \V 6 



N+ 1 1 _ 2tt{N + 1)L 

A Ar ~~ NA 
N+l 1 (N + 1)L 



A Alp Ay/1 - eL 2 



Note that the ratio m*/m is the same as in Eq. (|7.14J) for the gravitons. However, the 
mass terms are different and can be written symbolically in the form 



M = m 4 



12 

Xr Xr Xr 



xl o 



Xl 

xl 



■1 1 
•1 f 



+ im ■ 



12 

Xr Xr Xr ' " 

xj> 

o -l i 

xl -11 



In order to find the mass eigenvalues we apply a bi-unitary transformation meaning that 
left- and right-handed fields transform differently. Thus the states x an d the mass eigen- 
states ip are related by 



j 

XI 



Xr 



Xr, 



1 



N 



n 



n=l 



VN + T 

N-l 



N 



N 



sj N{N + 1) 
-^£exp(-2 7 rz. J ^ + 

v n=l 



(7.35) 



1 



VN + T 



V^(iV + 1 



Expressed by the mass eigenstates ip the action S for iV + 1 4D Dirac fermions reads 



S 



N 

3=0 

N-l 

- y2 I d 4 x ■ iplip R [m* + im(e~ 2m % - 1)] + ijjfyii™'* - im(e 2m % - 1)] 



iV 



From this result we read off that there is one massless fermion ip°, one heavy fermion ip 
with mass m±VN + 1 and iV — 1 fermions ip 1 , . . . ,ip N_1 with complex masses, whose 
squared absolute values read 



7m, 



m+ + im(e~ ^ — 1)| = m 4 + 4m sin (^y - ) + 2m*msin 



,2irn. 
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In contrast to the graviton case from Eq. ()7.29)1 . here we find an additional interference 
term 2m*msin(^p) in the mass spectrum. Since this looks a bit strange we propose a 
slightly modified discretisation procedure for the angular direction. Instead of d 6 x ~^ 
{x^ +1 ~ X'O/Ay 9 from Eq. (|7.34J) we choose the following prescription on a formal level 

d 6 -> % ■ - 



Ay? 

This does not change anything with the zero mode or the heavy mode, but the transfor- 
mations in Eqs. (|7.35J) lead to modified complex masses for the modes ip l . . . ip N ~ l . Their 
squared absolute values read in this case 

• . 7TTZ . . Q r> o o , 1XTX . 

|m* + zm(2sm — )\ =m^ + Am sin \~^)- 

Here we see that the interference terms from above do not exist and the fermion mass 
spectrum has exactly the same structure as the gravitons. 



7.7 Small Fermion Masses 

Here, we finish this chapter with a simple application, where the results for the fermions 
on the discretised curved disk can be applied directly to generate small fermion masses. 
For this purpose let us assume that our known 4D world is located on the centre site of 
the disk. At this place the particle standard model may couple to the 4D component Xr 
of the 6D Dirac field via a Yukawa coupling schematically given by L(H)xr, where L 
is a left-handed lepton doublet and (H) is the VEV of the Higgs doublet, respectively. 
In addition we consider a large number N of lattice sites so that the heaviest mode Q N 
decouples due to its large mass m^y/N + 1. From Eq. ()7.35|) it follows that the right- 
handed fermion Xr on the centre site essentially consists only of the zero-mode ip^ with 
a weight factor l/y/N + 1. Thus the Yukawa coupling of the left-handed SM doublet L 
with the right-handed fermion Xr leads to a suppression of the SM neutrino mass, 

L{H)xl - -^=MH)<- 

For example, the number iV ~ 10 24 might explain the lightness of the SM neutrinos [86 . 
For related applications in this context see also Ref. [S3]. Finally, we mention that this 
mass suppression mechanism can be considered as a discrete version of the wave function 
suppression mechanism in continuous higher dimensions from Ref. [HHj . 
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In this thesis we have investigated the dark energy problem from two different perspec- 
tives that both involve vacuum energy, or equivalently the CC, as the major dark energy 
candidate. In the first part we have started with scaling laws for the CC originating from 
RGEs in quantum field theory and quantum gravity. However, since these theories do not 
fix the physical interpretation of the corresponding renormalisation scale, some typical 
scales occurring in the context of cosmology have been chosen to fit into this role. In this 
sense we have considered the Hubble rate in addition to the scales characterised by the 
cosmological event and particle horizons. In the end, this leads to a time- dependent CC 
implying a non-trivial scaling of the matter energy density with the cosmic scale factor. 
Solving Einstein's equations therefore becomes more complicated, but the solutions we 
found exhibit some very interesting features. 

In the late-time epoch of the cosmological evolution we have found, apart from the well 
known de Sitter final states, also power-law and super-exponential solutions for the scale 
factor. Additionally, we have observed in some cases future singularities of the big rip and 
big crunch type, which usually appear only in more exotic dark energy models. In detail 
we have discussed in this analysis all combinations of scaling laws and renormalisation 
scale identifications and solved Einstein's equations analytically and numerically. As a 
result we have determined the fate of the universe for each case and studied the influence 
of the parameters on the solutions. By doing this it happened that some final states occur 
only for certain scaling laws or scale interpretations. Finally, the running of Newton's 
constant was taken into account in one case, which we treated numerically in more detail. 
In conclusion, the results found in this thesis feature several very different solutions for 
the cosmological late-time evolution. Based on this outcome the analysis may help to 
discriminate between different combinations of scaling laws and scales. 

The second part of the work consists of an investigation of vacuum energy in higher di- 
mensions, where we make use of the fact that the zero-point energy of quantum fields 
depends on boundary conditions. Starting from the Casimir effect in a continuous higher- 
dimensional space-time, where the Casimir energy density was determined and compared 
with results from the literature, we have adjusted the calculations to the case of discrete 
EDs. Within a detailed analysis we have discussed the influence of boundary conditions 
and the number of lattice sites on the value of the Casimir energy density, which con- 
tributes to the effective 4D CC. Special emphasis was placed on the effect of bulk field 
masses, which open up the possibility to considerably suppress the corresponding vacuum 
energy contribution. In addition, we have discussed the differences between the bosonic 
and fermionic case and observed some typical lattice artefacts for fermions. 

A nice motivation for discrete EDs comes from the model building sector in the form 
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of deconstruction. Representing its most important property, this framework describes 
discretised EDs on the basis of 4D quantum fields thereby avoiding problems that ap- 
pear in continuous higher dimensions. As quantum fields in this setup are not subject to 
boundary conditions their zero-point energy is a priori completely unconstrained. At this 
stage we have proposed a prescription to solve this problem by employing the correspon- 
dence between the discrete 5D setup and deconstruction. Therefore the vacuum energy 
of all 4D fields in deconstruction has been identified with the finite Casimir energy of the 
corresponding 5D quantum field. Within a specific model we have demonstrated this idea 
on a formal level for bosons and fermions. By using the suppression mechanism of large 
field masses in this setup, we have also shown that the vacuum energy contributions lie 
below or around the observed value of the CC. 

Even if the strong condition of tiny CC contributions might be satisfied for arbitrary small 
EDs by implementing sufficiently large bulk masses, we have found that this might not be 
true for discrete gravitational EDs. In these models there usually exists a strong coupling 
scale that depends on the size of the ED thereby implying an upper limit for the mass 
scales in the theory. Due to this UV/IR connection one has lost the ability to arbitrarily 
suppress the Casimir energy by considering large bulk masses. In this context, a lower 
limit on the ED size emerges. By applying our results for the Casimir effect we have 
found the minimal size R m i n to be roughly in the range (10 12 GeV) -1 . . . (10 7 GeV) -1 . As 
a consequence, one can exclude, e.g., Planck scale sized EDs in this framework. 

Motivated by the appearance of massive gravitons in discrete gravitational dimensions our 
last subject of investigation consists of a 6D model, where both EDs form a discretised 
two-dimensional curved disk. After discussing the disk in the continuum, we have applied 
a special discretisation procedure that features 4D gravitons with an interesting mass 
spectrum. The explicitly performed calculation of the masses has led, for instance, to the 
observation of a gap between the zero mode and the finite KK tower. As the main result 
of this subject, both mass scales appearing in this spectrum were found to be completely 
adjustable by the parameters of the curved disk. In contrast to a flat disk, here it is 
the curvature of the disk that makes this flexibility possible. In the future, it will be 
interesting to explore some applications and the strong coupling scale in this setup or 
more generalised scenarios [86J. As a foretaste, we have finally shown in this work the 
explicit implementation of fermions, which could be directly applied to generate small 
fermion masses. 

As our final conclusion let us mention that, despite many possible candidates for the source 
of the accelerated expansion of the universe, the CC should be the first one to look at. We 
have seen in this thesis that quantum effects lead not only to a dependence on boundary 
conditions as given by discrete EDs, but they may also induce a time- dependence. By 
investigating both subjects we have found new insights for cosmology and the structure 
of our space-time. 
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